NASA-CR-1 96093 


/ ' /' 6 

50 f> 


Monolith Catalysts for Closed-Cycle Carbon Dioxide Lasers 


NASA Grant NAG-1 -1051 


Final Report 

Grant Ending Date March 31, 1994 


submitted by 

Principal Investigator: 

Richard K. Herz 

Department of Applied Mechanics and Engineering Sciences 
University of California, San Diego 
La Jolla, California 92093-0310 

to 

Grant Technical Officer: 

David R. Schryer 
NASA Langley Research Center 
Hampton, Virginia 23665 


(NASA-CR— 196093) MONOLITH 
CATALYSTS FOR CLOStD-CYCLE CARBON 
OTOXIOE LASERS Final Report, period 
ending 31 Mar. 1994 (California 
Un i v. ) 50 p 


N94-35489 
Unc 1 as 


G3/36 0013756 



"Monolith Catalysts for Closed-Cycle 
Carbon Dioxide Lasers" 


TABLE OF CONTENTS 

ABSTRACT ii 

FINAL REPORT 1 

OBJECTIVE 1 

MANUSCRIPT 2 

APPENDIX I (program listing) 22 

APPENDIX II (program listing) 35 

(last page of final report) 47 



ABSTRACT 


The general subject area of the project involved the development of solid catalysts that have high 
activity at low temperature for the oxidation of gases such as CO. The original application 
considered was CO oxidation in closed-cycle CO2 lasers. The scope of the project was 
subsequently extended to include oxidation of gases in addition to CO and applications such as air 
purification and exhaust gas emission control. The primary objective of the final phase grant was 
to develop design criteria for the formulation of new low-temperature oxidation catalysts utilizing 
Monte Carlo simulations of reaction over NASA-developed catalysts. This work resulted in a 
paper published in the Journal of Catalysis . 
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FINAL REPORT 


OBJECTIVE 

The primary objective of the grant was to develop design criteria for the formulation of new low- 
temperature oxidation catalysts utilizing Monte Carlo simulations of reaction over NASA- 
developed catalysts. This work resulted in a paper published in the Journal of Catalysis in. A 
copy of this paper follows and serves as the main body of this final report. The computer 
programs used in this work are listed in two appendixes to this final report. 
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Two-Component Catalysts for Low-Temperature CO Oxidation: 

A Monte Carlo Study 

R. K. Herz,*- 1 A. Badlani,* D. R. Schryer,! and B. T. Upchurch! 

Chemical Engineering Group, Department of Applied Mechanics and Engineering Sciences, 
University of California . San Diego. California 92093-0310; and "tNASA Langley Research Center, 

Hampton , Virginia 23681-0001 

Received August 20, 1992; revised January 5, 1993 

The reaction of CO and 0> at low temperature over composite, noble-metal/reducible-oxide 
catalysis is simulated using Monte Carlo techniques. High activity for CO oxidation can be obtained 
over a composite material composed of a highly interspersed mixture of one type of site that 
adsorbs CO and 0 2 and another type of site that adsorbs 0 2 without significant CO inhibition. For 
example, the rale over Pd under \% of an atmosphere of CO at room temperature is predicted to 
increase 10 orders-of-m 2 gnitude with addition of \% of surface sites which adsorb 0 2 but not CO. 
For most reaction rules and parameter values, a roughly 50-50 mixture of the two types of sites 
gives the greatest activity per unit total surface area. This result is determined by the reaction 
stoichiometry and the fact that the two reactants primarily adsorb separately on the two different 
types of sites. In a randomly distributed mixture, the two types of sites have widely differing 
activities which depend on the local site configurations. The local site configurations of the most 
active sites in a random surface arc similar to site configurations found in a search for optimal 
configurations. The site configurations found in the search for optimal configurations were about 
20% more active than the random surfaces of the same overall composition. This relatively small 
increase may be due to the simple steric requirements of CO and 0 2 adsorption. We expect that 
similar searches for optimal site configurations will be more fruitful for reactions w ith more complex 
steric requirements, e iv*.i Academic prc*v Inc. 


INTRODUCTION 

For most of us who are familiar with con- 
ventional CO oxidation catalysts which are 
active only at temperatures above about 
I50°C, observation of rapid CO oxidation at 
room temperature can be startling. In our 
laboratories, for example, we have mea- 
sured, at room temperature, 56% conver- 
sion of a stoichiometrically balanced, atmo- 
spheric-pressure mixture of 1% CO and 
0.5% 0 2 in N 2 flowing at 1 .2 cnrVs over 0.2 g 
of Au/Mn0 2 (/). Au/Mn0 2 and related ma- 
terials are some of the few heterogeneous 
catalysts with activities at ambient condi- 
tions which may approach those of biologi- 
cal catalysts.' The objective of the work 
presented here is to help develop an under- 

1 To whom correspondence should be addressed. 
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standing of how these heterogeneous cata- 
lysts arc able to oxidize CO at ambient and 
subambient temperatures. 

There are several applications in which 
the catalytic oxidation of CO near ambient 
temperature is desirable. One such applica- 
tion is removal of CO as a contaminant of 
breathing air in enclosed spaces such as sub- 
marines and space vehicles and in burning 
structures or mines (2). Another is oxidation 
of CO in automobile engine exhaust during 
cold starts (2). A third is regeneration of 
C0 2 in transversely excited atmospheric 
pressure (TEA) C0 2 lasers ( 4 , 5). Stoichio- 
metrically balanced mixtures of CO and 0 2 
arc generated during the operation of C0 2 
TEA lasers through the decomposition of 
C0 2 by the electrical discharge that initiates 
the lasing process. Consumption of C0 2 and 
buildup of 0 2 degrades the performance of 
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such a laser. A CO oxidation catalyst that 
operates under low-temperature conditions 
is desirable for the development of closed- 
cycle or sealed C0 2 TEA lasers, which have 
applications such as mapping of earth’s 
wind patterns from space by laser Doppler 
anemometry and remote sensing of environ- 
mental pollutants by infrared spectrometry 
( 6 ). 

The reaction between CO and 0 : pro- 
ceeds at a negligible rate over metal cata- 
lysts near ambient temperature, even 
though the reaction is thermodynamically 
favored and there are several metals, includ- 
ing Pt and Pd, which can chemisorb both 
species and, thus, should permit a Lang- 
muir-Hinshclwood-type reaction mecha- 
nism to occur (7). The explanation fre- 
quently given for the lack of reaction under 
these superficially favorable conditions is 
that CO requires only a single vacant ad- 
sorption site and 0 2 requires two adjacent 
vacant sites, so that a clean surface exposed 
to a mixture of CO and 0 ; quickly becomes 
covered with CO, which prevents 0 2 ad- 
sorption, except at low C0-to-0 : ratios (8). 
Only at elevated temperatures, where CO 
desorption becomes significant, do enough 
adjacent adsorption-site pairs become avail- 
able to allow significant 0 2 chemisorption 
and reaction to occur. 

A demonstration of the foregoing expla- 
nation w'as presented by Ziff et al. (#)■ Using 
a simple stochastic (Monte Carlo) model 
which assumed equal sticking probabilities 
for CO and 0 2 and neglected CO desorption, 
they showed that chemisorption and reac- 
tion of both CO and 0 2 could occur on a 
uniform surface only when the ratio of the 
partial pressures of CO to 0 2 , P C o^o 2 < fell 
betw'een 0.59 and 1.0. At higher P co^o 2 > 
including the stoichiometrically balanced 
case where P C o lp o, = 2, the surface be- 
comes covered with CO. If the slicking 
probability of CO is greater than that of 0 2 , 
which would be the case for a Pt surface, 
these partial-pressure ratios would be 
shifted downward and further away from 
stoichiometry. 


The treatment by Ziff et al. is quite gen- 
eral and should apply to any single-compo- 
nent catalyst in which a species requiring 
only a single adsorption site is to be oxidized 
by 0 2 . A number of other workers have 
studied the CO + 0 2 reaction, or the general 
A + B, reaction, over surfaces using Monte 
Carlo (9-18) and cellular automaton simula- 
tions (19). Ertl and co-workers have used 
cellular automaton simulations to study the 
participation of Pt surface reconstruction in 
rate oscillations during CO oxidation (20, 
21). Related Monte Carlo studies have been 
made of the general A + B reaction over 
surfaces (22—30). Monte Carlo simulations 
are required in order to study the kinetics of 
surface reactions in which adsorbed species 
are not randomly distributed over the sur- 
face; traditional kinetic models using alge- 
braic rate equations assume such random 
distribution. 

The fact that a CO oxidation catalyst must 
perform tw-o distinct functions, adsorption 
of CO and dissociative adsorption of 0 2 , 
suggests that a “composite material” com- 
posed of an intimate mixture of two different 
materials, each with optimal properties for 
one of the tw'o catalytic functions required, 
would have significant low-temperature ac- 
tivity (_?/). In such an ideal “composite cata- 
lyst,” CO and 0 2 w'ould not compete for the 
same adsorption sites and, thus, CO inhibi- 
tion of low-temperature CO oxidation would 
be eliminated. 

A class of composite materials is cur- 
rently underdevelopment for application as 
low-temperature CO oxidation catalysts in 
COt TEA lasers and in breathing air purifi- 
cation. We have given the term “noble- 
metal/reducible-oxide” (NMRO) catalysts 
to this class of catalytic materials (32). Ex- 
amples include Au/Mn0 2 (33), Au/Fe 2 0 3 
(34), and Pd/Sn0 2 (35, 36). These materials 
have greater activity and stability than the 
conventional low-temperature CO oxida- 
tion catalyst, Hopcalite, which is a mixture 
of Cu and Mn oxides with small quantities 
of other oxides (37). In each of the NMRO 
catalysts, neither of the separate compo- 
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nents have appreciable activity for CO oxi- 
dation at low temperature and, thus, a syn- 
ergistic interaction is present in the two- 
component materials. Although the mecha- 
nisms of CO oxidation over these materials 
are not known in detail (32), a probable gen- 
eral explanation for their low-temperature 
activity is that the noble metal chemisorbs 
CO and the reducible oxide provides sites 
that dissociatively adsorb 0 2 . There are also 
indications that surface hydrogen or hy- 
droxyl participates in the dissociative ad- 
sorption of 0 2 and the oxidation of adsorbed 
CO at low temperature (J£). 

One barrier to our understanding of the 
mechanism of CO oxidation over these com- 
posite catalysts is that we do not know the 
surface structure and composition of these 
materials. Discrete noble-metal particles are 
present in many catalyst formulations. Re- 
cently, Brosilow el aL (9) adapted the model 
of Ziff el aL ( 8 ) to investigate the effect of 
a reducible-oxide support on the oxidation 
of CO over noble-metal particles. The effect 
of the reducible-oxide support was simu- 
lated by holding the coverage of adsorbed 
oxygen at saturation at the perimeter of no- 
ble metal particles. Near a stoichiometric 
ratio of CO and 0 : , PqqIPo , ~ 2. 1 , the reac- 
tion rate was nonzero for a distance of three 
or four noble-metal lattice sites with pene- 
tration from the perimeter into the noble- 
metal particle; the rate would be zero every- 
where on the particle in the absence of the 
reducible oxide at the perimeter. At lower 
CO to 0 2 ratios, the effect of the reducible 
oxide penetrates further into the noble- 
metal particle. These results suggest that 
very small noble-metal particles would pro- 
vide the highest rate per noblc-mctal surface 
atom, especially for stoichiometric mix- 
tures. 

Enhanced reaction at the perimeter of no- 
ble metal particles in contact with a reduc- 
ible-oxide support is not the only possible 
reason for the enhanced activity of NMRO 
catalysts. The activity of Pt/SnO : is in- 
creased by reducing pretrealments, and 
Hoflund and co-workers (39, 40) reported 


that such prelrealment leads to the forma- 
tion of Pt-O-Sn, Pt(OH) 2 , metallic Pt and 
Sn, and Pt-Sn alloy. In related work. Oh 
and Carpenter (41) found that a Pl-Rh cata- 
lyst which contained bimetallic Pt-Rh parti- 
cles had enhanced activity for CO oxidation. 
They proposed that Rh provided sites favor- 
able for 0 2 adsorption near CO adsorbed on 
Pt. Rh may have been present as a reducible 
oxide in the Pt-Rh bimetallic particles, mak- 
ing this catalyst a member of the NMRO 
class of materials. Logan and Paffett (42) 
studied the CO oxidation activity of a 50-50 
Pd-Sn surface alloy. They found that Pd-Sn 
had higher activity than pure Pd but also 
found that the Sn oxidized to form a partial 
SnO v overlayer. Work with Rh/Ti0 2 has 
demonstrated that the unique activity of this 
material after high temperature reduction 
pretreatment is related to the presence of 
TiO v overlayers covering much of the sur- 
face of the Rh particles (43, 44). Although 
the unique activity in these cases may result 
from reaction at the interface between the 
noble-metal and oxide overlayers, this in- 
ferfacc may be complex and the noble-metal 
and reducible-oxide components highly in- 
terspersed. 

In this work we investigate the reaction 
of stoichiometric mixtures of CO and 0 2 
over highly interspersed mixtures of noble- 
metal and reducible-oxide sites using Monte 
Carlo simulations. We demonstrate that a 
two-component catalyst, in which one com- 
ponent chemisorbs CO and the second com- 
ponent chemisorbs 0 2 but does not signifi- 
cantly chemisorb CO, makes possible the 
low-temperature reaction of stoichiometric 
mixtures of CO and 0 2 . We suggest that this 
feature accounts for the pronounced low- 
temperature activity of NMRO catalysts 
which readily catalyze the reaction of stoi- 
chiometric mixtures of CO and 0 2 at room 
temperature and below (4, 34). Monte Carlo 
simulations are performed in order to inves- 
tigate the effects of different 0 2 adsorption 
rules, the effect of changing the 0 2 adsorp- 
tion probability, the kinetic orders of the 
reaction, and the effects of changing the CO 
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desorption probability. Finally, a search is 
made for optimal configurations of surface 
sites as a first step towards the rational de- 
sign of composite catalysts. With the excep- 
tions of the study of uniform surfaces with 
heterogeneous boundaries by Brosilow ct 
al. (9) and the study of uniform surfaces 
with defects by Vlachos et al. {15 , 45), this 
is the first Monte Carlo study, of which we 
are aware, of catalytic reaction over nonuni- 
form, multicomponent surfaces. 

REACTION MODEL 

Noble-metal sites which adsorb CO are 
referred to in this work as a sites and the 
reducible-oxide sites which adsorb 0 2 are 
referred to as ft sites. Previous Monte Carlo 
simulation work on uniform surfaces has re- 
ferred to the CO + 0 2 reaction as a member 
of a general class represented by the reac- 
tion A + B 2 . The labels a and have been 
used previously to distinguish between dif- 
ferent sites on catalyst surfaces identified 
during thermal desorption experiments {46). 
Thus, the cases studied here can be referred 
to as examples of the general A + B : reac- 
tion over a two-component, a-fi surface. 

The identity of individual ft sites is not 
specified in the Monte Carlo simulations. In 
Pt/Sn0 2 , for example, actual individual (3 
sites might consist of one or more of the 
following: (Sn), , (SnH,)., (SnO^)., (SnO* 
OH v ) z , (SnOjfHj.);, (SnOH v ). , etc. Inter- 
spersion of noble-metal and reducible-oxide 
materials should lead to modification of the 
properties of each material. Such modifica- 
tions are not explicitly considered here, al- 
though they could be studied by determining 
the effect of varying the adsorption, desorp- 
tion and reaction probabilities assigned to 
each site. Rather, these simulations mainly 
probe geometric effects such as the effects 
of changing the relative positions of the two 
sites in a surface. 

The main loop in the simulation program 
is shown in Fig. 1. This program was devel- 
oped in order to consider sets of parameter 
values where adsorption and reaction prob- 
ability values are within several orders-of- 


magnitude of each other. The adsorption, 
desorption, and reaction probabilities are in- 
dependent of surface coverage. The 0 2 
sticking probability is the same for all al- 
lowed site pairs, and the reaction probability 
is the same whether the O atom is adsorbed 
on an a site or a /3 site. Surface diffusion is 
not described. A lattice of square sites was 
used in all simulations. 

The following “base set” of parameter 
values was used below, except where speci- 
fied otherwise: 

(a) CO sticking probability (probability 
that a gas phase CO molecule striking a va- 
cant a site will adsorb): p { = 1; 

(b) 0 2 sticking probability (probability 
that a gas phase 0 2 molecule striking a va- 
cant pair of sites that correspond to the spec- 
ified set of 0 2 adsorption rules — e.g., a-a, 
a- j3 will adsorb): p 2 = 1; 

(c) CO desorption probability (probabil- 
ity that a CO molecule will desorb when a 
site is selected randomly at a frequency 
equal to the CO collision frequency and the 
site is occupied by an adsorbed CO; as de- 
fined, this probability is inversely propor- 
tional to CO pressure): p } = 0; 

(d) reaction probability (probability that 
a reaction event will occur when a pair of 
neighboring sites is selected randomly and 
the site pair is occupied by one adsorbed 
CO molecule and one adsorbed O atom): 
Pa = 0.4; 

(e) actions allowed on diagonally adja- 
cent pairs of sites, 

(f) stoichiometrically balanced ratio of 
gas-phase CO and 0 2 : P c q/Po 2 = 2- 

The differing molecular weights of CO and 
0 2 were taken into account when determin- 
ing the relative collision frequencies of these 
reactants in the simulations. 

Reaction rates are reported relative to the 
CO collision frequency. The results pre- 
sented here are expected to correspond 
most closely to reaction over Pt/Sn0 2 and 
Pd/Sn0 2 catalysts at room temperature, 
where the CO sticking probability is near 
one and the CO desorption probability is 
near zero. The results presented also apply 
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Fig. 1 . Main loop in the Monte Carlo simulations, r = new random number between 0 and l 
generated at each step within the loop. V| = fraction of iterations that are checks for reaction — >s 
(^co.basc^coVH + >*2 {P a where {P C0 ^ Ast / P co ) is the ratio of the CO partial pressure for the 
base set of conditions to the CO partial pressure for the specific simulation run. v, is defined such 
that the frequency of checks for reaction in the program at the base set of conditions equals the 
frequency of CO collisions and such that the frequency of checks for reaction in “real time" is 
independent of CO pressure. v 2 = fraction of gas molecules colliding with the surface that are CO 
molecules, p, = CO sticking probability. p : - 0 2 sticking probability. p 3 = CO desorption probability. 
Pa = reaction probability. See text for further explanation of probabilities. 


to CO sticking probabilities less than one 
when the other probabilities and the reac- 
tion rate are normalized properly with re- 
spect to the CO sticking probability. Here, 
we refer to the unnormalized probabilities 
and rates for simplicity. 

All simulations were performed starting 
from a clean surface and run to steady-state 
conditions. In most cases, steady state was 
reached within 1000 Monte Carlo steps, 
where, in one Monte Carlo step, the loop 
in Fig. 1 is repeated the number of times 
equal to the total number of sites in the 
array. Cases which did not reach steady 
state within about 1000 Monte Carlo steps 
eventually deactivated completely. AJ1 reac- 
tion rates reported are average steady-state 
rates. In most cases, these rates were ob- 
tained by starting averaging after attainment 
of steady state at 3600 Monte Carlo steps 
into the run and averaging over the subse- 
quent 3900 Monte Carlo steps. 


For the parameter values used here, the 
effect of /3 sites on a sites is localized and 
propagates only a short distance into a 
"patch” or "particle” of a sites. Figure 2 
shows the variation in reaction rate on a 



Column numbef ol a sHes Irom the a- 3 interface 

Fig. 2. Local rate vs distance from a semi-infinite 
linear a-/3 interface. The local rate is equal to the num- 
ber of reaction events on an a site per CO collision 
with that site. 
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sites with distance from a semi-infinite a-(3 
interface for the base set of parameters and 
for a case in which 0 2 can adsorb on any 
pair of neighboring sites. The reaction rate 
is zero everywhere except for the three col- 
umns or rows of a sites nearest the a-/3 
interface. This result is similar to that ob- 
tained by Brosilovv et al. (9) at P c 0 /P 0l = 
2.1 for adsorption-limited reaction. 

RANDOM a-p SITE DISTRIBUTIONS 

A 42 x 26 lattice with 1092 total sites 
and periodic boundaries was used for the 
simulations of all-a and random a-/3 sur- 
faces. Random a-p surfaces represent the 
opposite extreme in geometric a-/3 config- 
uration from the noble-metal particles on 
reducible-oxide supports simulated by 
Brosilow et al. (9). Some Pt-Sn surface 
alloys have a highly interspersed geometry 
(42); reaction over other Pt-Sn alloys may 
involve an intermediate a-j 3 geometry, 
which might be obtained in simulations by 
“annealing'' initially random surfaces. 

For the case in which 0 2 can adsorb on 
any site pair and the CO desorption proba- 
bility is zero, the reaction rate for the all-a 
surface is zero. The rate is nonzero for the 
all-a surface for nonzero CO desorption 
probabilities. For low CO desorption proba- 
bilities, the rate increases dramatically 
when p sites are added to the surface. The 
relative increase is inversely related to the 
CO desorption probability. At ambient tem- 
perature over Pd, the CO desorption fre- 
quency from a site is on the order of 10° s" 1 
(7). This desorption frequency leads to a CO 
desorption probability, as defined above, of 
10' 7 at P Q o = 0.01 atm, a CO pressure char- 
acteristic of C0 2 TEA lasers (4). With this 
value of the CO desorption probability, the 
overall rate increases by 10 orders-of-mag- 
nitude — from 10“ M reaction events per CO 
collision with the surface to 10~ 4 — as the 
fraction of P sites is increased from 0 to 
1%. This result demonstrates the extreme 
sensitivity of steady-state CO oxidation 
rates over Pd and Pt to the addition of small 
amounts of an oxygen adsorbing component 


that is not inhibited by CO. In an experimen- 
tal demonstration of a related phenomenon, 
Mundshau and Rausenberger (47) used pho- 
toelectron microscopy to show the presence 
of defect sites present in low concentration 
in the surface of single crystal Pt that adsorb 
CO only weakly and initiate the transient 
burn-off of inhibiting CO overlayers by 0 2 . 

A CO desorption probability of zero was 
used in most of this work in order to simulate 
the strong CO adsorption on noble metal 
sites at ambient temperatures and, thus, the 
strong inhibition of 0 2 adsorption on Pt- 
group noble metals. In general, the absence 
of reactant desorption and surface diffusion 
in Monte Carlo simulations produces results 
which highlight the geometric effects pres- 
ent in the physical system modeled. 

The reaction rate is zero for all a-plus-/3 
surfaces with no CO desorption when 0 2 
adsorption is not allowed on /3-/3 pairs. This 
is because all a sites are saturated with CO 
in this case, thus blocking 0 2 adsorption on 
a-a and a-fi pairs. The rate is nonzero with 
no CO desorption on a-plus-/3 surfaces 
when oxygen adsorption is allowed on /3-/3 
pairs. 

We focus on the two cases of adsorption 
of 0 2 on (a) any site pair and (b) only fi-fi 
site pairs. In each of these two cases, the 
0 2 sticking probability is the same on all 
aliowed site pairs. In real catalysts, we 
would expect that the sticking probability 
would be different on different site pairs. 
The two cases considered here are limiting 
cases, with real systems possibly having be- 
havior intermediate between these two 
cases. 

Figure 3 shows the steady-state rate-per- 
a-sile vs the fraction of p sites present in 
the surface (“fraction-/}” below). The rate- 
per-a-site is equal to the average number of 
reaction events per CO collision with an a 
site. In many cases, the noble metal compo- 
nent represented by the a sites may be the 
major cost factor in a catalyst. The open 
circles are for the case in which 0 2 can ad- 
sorb on any site pair. The solid diamonds 
are for the case which 0 2 can only adsorb 
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Fracl>on of 3 s?es 

Fig. 3. Steady-slate rate* per-a- site vs fraction-/? for 
the base set of parameters. The rate-per-a-sitc is equal 
to the average number of reaction events per CO colli- 
sion with an a site. 


on (3-13 pairs. At each nominal fraction-/?, 
20 runs were performed to show the varia- 
tion in rate from run to run, except for only 
three runs at fraction-/? = 0.99 for the solid 
diamonds. At each specific fraction-/?, the 
variation in rate is caused by the fact that 
different random configurations of a and (3 
sites have somewhat different activities. 

For the case in which 0 : adsorption is 
allowed on (3~{3 pairs only, the rate-per-a- 
site is zero at fraction-/? = 0. since no CK 
can adsorb, and increases continuously as 
the fraction of (3 sites is increased. The rate- 
per-a-site is somewhat more linear than sim- 
ply being proportional to the square of the 
fraction of (3 sites. The rate-per-a-site ap- 
proaches a value of 0.33 reaction events per 
CO collision with an a site as the fraction 
of (3 sites approaches one, that is, in the 
limit of isolated a sites surrounded by (3 
sites. In surfaces with high fraction-/? but 
before this limit, two widely separated a 
sites are more active than two adjoining a 
sites but are less active than two a sites 
separated by a distance of one or two (3 
sites. In the latter case, 0 2 adsorption and 
reaction are enhanced in the region between 
the two a sites since there is an increased 
rate of formation of vacant /?-/? site pairs 
in this region. 

For the case in which 0 2 adsorption is 


allowed on any site pair, the rate-per-a-site 
goes through a maximum at fraction-/? = 
0.5. To the left of the rale maximum, the 
rate increases as more (3 sites are added to 
the surface because 0 2 can adsorb on (3 sites 
next to CO molecules adsorbed on a sites. 
To the right of the rate maximum, the rate- 
per-a-site decreases as more (3 sites arc 
added to the surface. This occurs because 
a sites are becoming increasingly dispersed 
among /? sites and, thus, are becoming more 
susceptible to deactivation by irreversibly 
adsorbed oxygen atoms. These oxygen 
atoms cannot desorb because the oxygen 
desorption probability is zero. They cannot 
be removed by reaction because there are 
no neighboring a sites in the local vicinity 
which are not also deactivated by oxygen. 

Figure 4 shows a “snapshot” of the sur- 
face for the two cases at fraction-/? = 0.5. 
The reaction rates are approximately equal 
for both cases at this fraction-/?. Note that 
there is less CO present on the surface for 
the case where 0 2 can adsorb on any site 
pair. This is because some a sites are cov- 
ered with adsorbed oxygen atoms. Rela- 
tively isolated a sites are deactivated by this 
adsorbed oxygen. Oxygen adsorbed on a 
sites near other a sites covered with CO can 
be removed by reaction. At fraction-/? = 
0.5, the deactivation of isolated a sites by 
oxygen roughly balances the contribution to 
reaction by oxygen adsorption on a sites 
which remain active. Qualitatively, ad- 
sorbed CO and O are interspersed for the 
case where 0 2 can adsorb only on (3-(3 pairs, 
whereas adsorbed CO and O are found pri- 
marily in separate patches for the case 
where 0 2 can adsorb on any site pair as a 
result of the coverage and deactivation of 
isolated a sites by adsorbed O. 

For the case in which 0 2 adsorption is 
allowed only on neighboring /?-/? pairs, iso- 
lated (3 sites are inactive, of course. Isolated 
(3 sites are also inactive for the case in which 
0 2 adsorption is allowed on any site pair 
when the CO desorption probability is zero. 
With a nonzero CO desorption probability 
for this second case, an isolated (3 site is 
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Vacant a site 

CO on a site 
Oxygen on a site 
Vacant p site 

Oxygen on p site 


Oxygen adsorption on p-p Oxygen adsorption on any pair 
pairs only 


Fig. 4. Instantaneous surface configurations during steady-state reaction for two cases of O; adsorp- 
tion. The fraction of £ sites is 0.5 for both cases and the base set of parameters is used. Only parts 
of the complete surfaces are shown. 


active but has a much lower activity than a 
pair of neighboring (3 sites. An isolated pair 
of / ? sites is active for both cases (0 : adsorp- 
tion only on /?-*/? pairs, 0 2 adsorption on 
any pair) with no CO desorption. Such an 
isolated (3-/3 pair, however, is more active 
for the case in which 0 : can adsorb on any 
site pair since 0 2 adsorption can involve a 
neighboring a site and one of the /? sites, 
not just the (3-(3 pair itself. The height of 
the solid curve above the dashed curve at 
low fraction-/? shows the extent to which 0 2 
adsorption on a-a and a~(3 pairs contrib- 
utes to the reaction rate for the case in which 
0 2 adsorption can occur on any site pair. 
The differences between these two cases 
decrease as the oxygen sticking probability 
decreases. 

For the case in which 0 2 adsorption is 
allowed on any site pair, note in Fig. 3 that 
the scatter of rates at a given fraction-/? is 
higher at large fraction-/? than at low 
fraction-/?. This scatter at large fraction-/? is 
caused by the fact that different proportions 
of site configurations susceptible to oxygen 
poisoning are generated at the start of the 
different runs. The surface becomes com- 
pletely oxygen poisoned for a fraction of 
(3 sites equal to 0.8 and above. The time 


required for complete deactivation to occur 
in these cases is much longer than the nor- 
mal start-up transient for runs in which the 
surface remains active. 

For the runs shown in Fig. 3, actions in- 
volving diagonally adjacent sites are al- 
lowed. When "diagonal actions"' involving 
these "diagonal pairs" are not allowed, the 
reaction rate decreases by about one-third. 
However, when diagonal actions arc not al- 
lowed there is no qualitative change in be- 
havior, with one exception: for the case in 
which 0 2 can only adsorb on /?-/? pairs, 
the rate-per-a-site drops to zero as the frac- 
tion of (3 sites approaches one rather than 
approaching a nonzero value. That is, indi- 
vidual a sites surrounded by / ? sites are inac- 
tive at steady state when diagonal actions 
are not allowed. After an initial transient in 
which CO adsorbs on these sites and reacts 
with oxygen atoms which have adsorbed on 
neighboring (3 sites, the reaction ceases. 
This is because oxygen is removed by reac- 
tion from the (3 sites to the left and right and 
top and bottom of an a site and this oxygen 
cannot be replaced. Since, at steady slate, 
the /? sites neighboring the a site diagonally, 
and essentially all other /? sites, are filled 
with oxygen atoms, these vacant /? sites to 
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Fig. 5. Instantaneous surface configuration during 
steady-slate reaction fora case in which 0 : can adsorb 
on pairs only; the fraction of 0 sites on the entire 
surface is 0.99, and the base set of parameters is used 
except that actions involving diagonally adjacent pairs 
of sites are not allowed. AH a sites shown arc covered 
by adsorbed CO, as indicated by the large black circles. 
The a site located in the upper left and the two a sites 
located on the far right are inactive because 0 : cannot 
adsorb on the vacant /3 sites (small black dots) located 
to their top and bottom and left and right. The two a 
sites located just left of center are active since 0 : can 
adsorb and react on the pair of £ sites directly between 
them. 


the left and right and top and bottom do not 
have any vacant adajeent 0 sites with which 
to form a vacant 0-0 pair to allow further 0 2 
adsorption. This inactive steady-stale site 
configuration is shown in Fig. 5. 

Figure 6 shows the same runs plotted in 
Fig. 3 but now plotted as the “overall rate" 
vs fraction-0. The overall rate is defined as 
the average number of reaction events per 
CO collision with any site on the surface. 
This overall rate corresponds to the rate that 
would be measured in the laboratory in 
moles per time per unit BET surface area. 
Note that both cases show a maximum in 
overall rate vs fraction-0. Although the rate- 
per-a-site is highest at high fraction-0 for 
the case in which 0 2 adsorbs only on 0-0 
pairs, the overall rate is relatively low be- 
cause of the low fraction of a sites. For both 
cases at low fraction-0 (0 to 0.3), the overall 
reaction rate in each case is roughly propor- 
tional to the fraction of a sites times the 
square of the fraction of 0 sites. The overall 
rate for the case in which 0 2 adsorbs only 
on 0-0 pairs continues to be roughly de- 


scribed by this proportionality from low to 
high fraction-0. 

The results shown in Fig. 6 arc in qualita- 
tive agreement with the experimental results 
of Upchurch et al. (48). They found that the 
CO oxidation activity, per unit weight, of a 
series of Pt/Sn0 2 catalysts was highest at 
intermediate Pt-Sn0 2 weight ratios. Sur- 
face composition measurements, especially 
at low noble-metal surface fractions, will be 
required in future studies in order to distin- 
guish between different 0 2 adsorption rules. 

Alternate Rules for 0 2 Adsorption 

Two different sets of rules for 0 2 adsorp- 
tion were discussed above: 0 2 adsorption 
on any site pair and 0 2 adsorption only on 
0-0 pairs. The comparison between these 
two sets of rules demonstrates that, when- 
ever 0 2 can adsorb on 0-0 pairs, the pre- 
dominate mode of reaction involves 0 2 ad- 
sorption on these pairs, with reaction 
involving 0 2 adsorption on a-0 and, espe- 
cially, a-a pairs contributing only slightly. 
The main contribution seen from allowing 
0 2 adsorption to involve a sites was to allow 
for oxygen poisoning of highly dispersed a 
sites at high fraction-0. 

We have also investigated the behavior 



Fig. 6. The same runs in Fig. 3, now plotted as the 
overall rate vs fraclion-/3. The overall rate shown here 
is equal to the average number of reaction events per 
CO collision with any site on the surface. This rale 
would be proportional to the observed rate per total, 
or BET, surface area measured experimentally over 
an actual catalyst. 


10 



228 


HERZ ET AL. 


with other sets of site pairs on which 0 2 
adsorption is allowed. The case in which 0 2 
adsorption is allowed on a-/? and /?-/? pairs 
behaves similarly to the case in which 0 2 
adsorption is allowed on all site pairs. 

Two sets of rules were investigated which 
do not allow 0 2 adsorption on /?-/? site pairs: 
0 2 adsorption only on (a) a-/? pairs and (b) 
a-a and a-/? pairs. For these sets of rules, 
the surfaces are inactive when the CO de- 
sorption probability is zero, since all a sites 
are covered by CO, preventing 0 2 adsorp- 
tion. As the CO desorption probability is 
increased slightly (e.g., to 0.01) from zero, 
the reaction rate becomes nonzero at most 
fractions of / ? sites. With a nonzero CO de- 
sorption probability, the ratc-pcr-a-site 
goes through a maximum at intermediate 
fraction-/?. The rate remains low at low 
fraction-/3 because of CO poisoning of the 
a sites in a-a and a-/? site pairs. In many 
cases, the rate remains zero at high 
fraction-/? because of oxygen poisoning of 
a sites. Specifically, for 0 2 adsorption on 
a-fi i pairs only and for a CO desorption 
probability of 0.01, the rate remains zero at 
fraction-/? = 0.7 and 0.9 for diagonal actions 
not allowed but is nonzero at these 
fraction-/? when diagonal actions are al- 
lowed. For 0 2 adsorption on a-a and a-/? 
pairs only and for a CO desorption probabil- 
ity of 0.01, the rate is zero at fraction-/? 
= 0.9 whether or not diagonal actions are 
allowed. The effect of allowing diagonal ac- 
tions is to enhance the reaction rate some- 
what, except in the case mentioned above 
for 0 2 adsorption on a-/3 pairs only where 
inactive surfaces become active if diagonal 
actions are allowed. 

Variation of 0 2 Sticking Probability 

In the simulations presented above, the 
0 2 sticking probability was set equal to the 
CO sticking probability, as in the Monte 
Carlo models of Ziff and co-workers (S, 49). 
Over noble metals, 0 2 sticking probabilities 
are lower than CO sticking probabilities 
(50). Here, the effect of reducing the 0 2 



Fraction ol p sites 

Fig. 7. Variation of 0 2 slicking probability, p : . 0 2 
can adsorb on any site pair and the base set of parame- 
ters was used except for the variation in p 2 . Note that, 
as p 2 decreases, the shape of the curve for this case 
approaches the shape of the curve in Fig. 3 for the case 
in which Oj can adsorb on pairs only. 

sticking probability is investigated for cases 
in which diagonal actions are allowed. 

For the case in which 0 2 adsorption is 
allowed only on /?-/? site pairs, decreasing 
the 0 2 sticking probability decreases the 
rate at all fraction-/?. For the case in which 
0 2 adsorption is allowed on any site pair, 
decreasing the 0 2 sticking probability de- 
creases the rate at low fraction-/? because 
of the reduced rate of 0 2 adsorption on all 
site pairs. This behavior is shown in Fig. 7. 
Decreasing the 0 2 sticking probability in- 
creases the rate at high fraction-/? because 
oxygen poisoning of highly dispersed a sites 
is reduced. The surface at fraction /? = 0.8 
goes from inactive to active when the 0 2 
sticking probability is decreased from 1 .0 to 
0.5. As a result of this behavior, the rale 
maximum shifts to higher fraction-/? and re- 
duces in amplitude as the 0 2 sticking proba- 
bility is decreased. The overall result is that, 
as the 0 2 sticking probability is reduced to 
lower values, the case in which 0 2 adsorp- 
tion is allowed on any site pair exhibits be- 
havior more similar to the case in which 0 2 
adsorption is allowed only on /3-/? pairs. 

Apparent Orders of Reaction 

The dependence of the reaction rate on 
changes in reactant pressure over small 
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ranges in reactant pressure can be reported 
in terms of a power-law rate expression. 
When the CO and 0 2 pressures are varied 
independently, the rate can be expressed as 
r = kP^oP o 2 . This rate expression is an em- 
pirical correlation of results and does not 
represent a kinetic mechanism, and the re- 
action orders are not restricted to integer 
values. For the Monte Carlo results pre- 
sented here, the apparent order of the reac- 
tion with respect to CO was determined by 
doubling the CO pressure, keeping the 0 2 
pressure constant and then using the 
equation 


0 = 


ln(r : V } ) 
In 2 ' 


0 ) 


where </> is the apparent order with respect 
to CO, r, is the rate (number of reaction 
events per unit “real time") at the base CO 
and 0 2 pressures in a stoichiometrically bal- 
anced ratio, and r 2 is the rate obtained when 
the CO pressure was doubled. The apparent 
order with respect to 0 : , y, was determined 
in a similar manner by doubling the 0 2 pres- 
sure and holding the CO pressure constant. 

For a stoichiometrically balanced mixture 
of CO and 0 2 , the rate can be expressed as 
r = kP £ o = k'PQ^ where P co = 0 .5P 0 , and 
where 77 is the apparent overall order. For 
the Monte Carlo results presented here, the 
apparent overall order of the reaction was 
determined by doubling both the CO and 0 2 
pressures while keeping the reactant pres- 
sure ratio stoichiometrically balanced. The 
overall order was then determined from the 
equation 


V = 


In (r : /r,) 
In 2 ’ 


( 2 ) 


where r, is the rate at the base CO and 0 2 
pressures in a stoichiometrically balanced 
mixture, and r 2 is the rate obtained when 
both the CO and 0 2 pressures are doubled. 
The overall order for a stoichiometric mix- 
ture will equal the sum of the individual 
orders, 77 = <f> + y, when the overall and 


individual orders are determined for infini- 
tesimally small variations in CO and 0 2 
pressures, or when the true kinetics obey 
the power-law expression exactly. 

There are only a limited number of experi- 
mental studies that report the kinetics of the 
CO oxidation reaction at low temperature 
over NMRO catalysts. Although these stud- 
ies have not provided a clear picture of the 
dependence of rate on reactant pressure, 
all show that the kinetics are different than 
those observed over single-component no- 
ble-metal catalysts. Stark and Harris (4) 
determined that the overall order of the 
reaction was approximately one for stoi- 
chiometrically balanced mixtures over Pd/ 
Sn0 2 and Pt/Sn0 2 , and Badlani (/) obtained 
the same result over Au/Mn0 2 . Over Pd/ 
Sn0 2 at low temperature. Bond ct aL (J5) 
found that the reaction was slightly negative 
order in CO. They also found that the reac- 
tion was approximately half order in 0 2 over 
a range of 0 2 concentration which depended 
on Pd concentration and temperature and 
tended to zero order at high 0 2 concentra- 
tions. Sampson and Gudde {51) studied stoi- 
chiometric mixtures of CO and 0 2 over a 
“precious metar'-Sn0 2 catalyst at low tem- 
perature and found that the reaction was 
zero order in CO and first order in 0 2 . Logan 
and Paffett {42) determined that the reaction 
was slightly positive order in both CO and 
0 2 over a 50-50 Pd-Sn surface alloy. 

Figure 8 shows the apparent reaction or- 
ders vs fraction-73 for the Monte Carlo simu- 
lation case in which 0 2 can adsorb on any 
site pair and for the base set of parameters. 
One interesting observation is that the over- 
all order for the stoichiometric mixture is 
roughly constant and equal to 0.5 over most 
of the range of surface composition. Al- 
though the separate orders with respect to 
CO and 0 2 change substantially over this 
range, these changes compensate each 
other. At small fraction-/?, reaction orders 
are similar to those that would be observed 
over a noble metal: the rate is positive order 
in 0 2 and negative order in CO due to CO 
inhibition of 0 2 adsorption on a sites. At 
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Fig. 8. Reaction orders vs fraction-^ for the case in 
which 0 2 can adsorb on any site pair and for the base 
set of parameters. 

large fraction-/?, the apparent orders are de- 
termined primarily by changes in oxygen 
poisoning of relatively isolated a sites: as 
the 0 2 pressure increases the poisoning be- 
comes more severe and the rate decreases; 
as the CO pressure increases a sites become 
more resistant to oxygen poisoning and the 
rate increases. 

Figure 9 shows how the overall order var- 
ies with the 0 2 sticking probability. When 
the 0 2 sticking probability is 0.01, the over- 
all order is approximately one, in agreement 
with most experimental measurements of 
reaction of stoichiometric mixtures over 
NMRO catalysts at low temperature (/, 4, 



Racoon of 0 ales 

Fig. 9. Effect of varying the 0 2 sticking probability, 
p 2 t on the overall order for the case in which 0 2 
can adsorb on any site pair and for base set of param- 
eters. 



Fig. 10. Instantaneous surface configuration during 
steady-state reaction for the case in which 0 2 can ad- 
sorb on any site pair; the fraction of /3 sites on the 
entire surface is 0.9, and the base set of parameters is 
used except that the 0 2 adsorption probability, p 2 , 
equals 0.01. 


51), and the reaction is first order in 0 2 
and zero order in CO, in agreement with 
the experiments of Sampson and Gudde 
(51). Figure 10 shows a characteristic con- 
figuration of a surface with 90% /? sites 
and an 0 2 sticking probability of 0.01. At 
any given time during steady state, the a 
sites are nearly saturated with CO mole- 
cules, /? sites neighboring a sites are nearly 
vacant, /? sites far from a sites are satu- 
rated with oxygen. These results predict 
that the reaction rate is limited by the 
adsorption of 0 2 at (or reoxidation of) 
reducible oxide sites located at the inter- 
face between the noble metal and reducible 
oxide components. 

Figure 11 shows the overall order vs 
fraction-/? for the case in which 0 2 can ad- 
sorb only on >3— >3 pairs and for the base set 
of parameters. Since there is no CO inhibi- 
tion or oxygen poisoning in this case, the 
order in CO decreases toward zero and the 
order in 0 2 and the overall order increase 
toward one as the 0 2 sticking probability 
decreases. At low 0 2 sticking probability, 
the behavior for the two cases of 0 2 adsorp- 
tion rules is similar since 0 2 adsorption on 
a-a and a-/? pairs becomes insignificant for 
the case in which 0 2 can adsorb on any site 
pair. 
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Frackjr of p sites 

Fig. 1 1 . Effect of varying the 0 2 sticking probability, 
/;>, on the overall order for the case in which 0 2 can 
adsorb only on fi-p pairs and for the base set of param- 
eters. 


Variation of CO Desorption Probability 

Oyer metal catalysts where CO inhibition 
dominates at low temperature, CO desorp- 
tion is a critical kinetic step in determining 
the rate. CO desorption is also the most 
highly activated step in the mechanism of 
CO oxidation over metals (7). At relatively 
high P CO IP 0 , , where the reaction is negative 
order in CO pressure, the overall reaction 
rate increases with temperature primarily as 
a result of an increased rate of CO desorp- 
tion and a lower inhibiting CO coverage. 
Thus, apparent activation energies of the 
overall reaction {42) are in the same range, 
80-120 kJ/mol, as the activation energy for 
CO desorption (7). 

CO desorption is less critical over NMRO 
catalysts than over metals because CO inhi- 
bition is less important. Over these NMRO 
catalysts at relatively high Pq 0 !P 0 , CO oxi- 
dation is only slightly positive {42) or 
slightly negative order in CO (J5) and the 
apparent activation energies are substan- 
tially lower — roughly 20-40 kJ/mol (/, 
42 ) — than for metal catalysts. 

The effect of increasing the CO desorp- 
tion probability from zero was investigated 
with the Monte Carlo model for cases in 
which diagonal actions are allowed. For the 
case in which 0 2 adsorption is allowed only 
on /?-/? site pairs, changing the CO de- 


sorption probability from zero to 0.01 pro- 
duced only a slight decrease in rate at all 
fraction-/?. Increasing the CO desorption 
probability further simply results in further 
decreases in reaction rate at all fraction-/?. 

For the case in which 0 2 adsorption is 
allowed on any site pair, changing the CO 
desorption probability slightly from zero to 
0.01 causes the all-a surface to become ac- 
tive but results in only slight changes in rate 
at higher fraction-/?. As the CO desorption 
probability increases, the rate at low 
fraction-/? increases and the rate at high 
fraction-/? decreases such that the rate maxi- 
mum shifts to lower fraction-/?. This trend 
of a shift of the rate maximum to lower 
fraction-/? continues as the CO desorption 
probability is further increased, as shown in 
Fig. 12. 

Continuing with the cases shown in Fig. 
12, for a CO desorption probability of 0.3, 
the surface at fraction-/? = 0.7 has become 
inactive as more a sites become susceptible 
to oxygen poisoning as the CO coverage on 
these sites decreases with the increase in 
CO desorption rate. Except for cases such 
as this in which a surface becomes com- 
pletely and irreversibly deactivated, the 
change in rate over a given surface with 
increase in CO desorption probability is re- 
versible when theCO desorption probability 
is later decreased. 



Fraction ot p sries 

Fig. 12. Variation of CO desorption probability, p y . 
0 2 can adsorb on any site pair and the base set of 
parameters was used except for the variation in p y . 
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SEARCH FOR OPTIMAL SITE 
CONFIGURATIONS 

One question that arises is what the opti- 
mal configuration of a and /} sites on the 
surface is in order to obtain the highest rate- 
per-a-site or the highest overall rate. One 
goal, after the correct mechanism is deter- 
mined, would be the ability to design and 
prepare the optimal catalyst. Because of the 
low-temperature operating conditions asso- 
ciated with many of the applications of 
NMRO CO oxidation catalysts, an optimal 
but thermodynamically unstable site con- 
figuration may be kinetically stable over op- 
erating periods of practical length. 

One possible way to search for the opti- 
mal surface configuration for a given set of 
rules and parameters would be to lake an 
active site configuration found from a run 
with a random surface and propagate it peri- 
odically across a larger surface. During each 
run the cumulative number of reaction 
events occurring on each site was recorded 
after the start-up transient. This allowed us 
to look at each surface configuration, search 
for the most active a site, for example, and 
then examine the site configuration sur- 
rounding this site. We found that the most 
active site configurations changed as the 
rules for 0 2 adsorption changed. Figure 13 
shows several most-active site configura- 
tions at various fraction -/} for cases in which 
0 2 can adsorb on any site pair and diagonal 
actions are allowed. Similar configurations 
and rates are obtained when diagonal ac- 
tions are not allowed. For each configura- 
tion, the most active a site is the central site 
in the surrounding 7x7 site array shown. 
Note that the local a-/} ratio and the rate 
on the most active a site are roughly con- 
stant even though the average a-/} ratio and 
the average rate change substantially be- 
tween different sets of patterns. Three of 
the five patterns shown have a local 
fraction-/? in the 3 x 3 array of sites centered 
on the most active a site = 0.44 (4/9), one 
pattern has a local fraction-/} = 0.56 (5/9), 
and one has a local fraction-/} = 0.33 (3/9). 

Another approach to finding an optimal 


site configuration would involve “evolu- 
tion” of a surface toward a more active con- 
figuration. Starting with a random surface 
of a desired fraction-/}, for example, the 
computer program could switch one of the 
least active a sites to a /} site and one of 
the least active /} sites to an a site. If this 
“mutation” resulted in an increase in rate, 
it would be preserved and another mutation 
would be tried. If the mutation resulted in 
a decrease in rate, it would be reversed. 

The approach we have pursued most ex- 
tensively involves searching for optimal pe- 
riodic configurations of a and /} sites. First 
we specify the size of a square “base 
array.” Second we choose a fixed number 
of a and /} sites that will populate this base 
array. Third, the specified number of a and 
/} sites are placed in the base array in one 
of the statistical combinations that are possi- 
ble. Finally, the base array is propagated 
periodically and equally in two dimensions 
for a specified number of repetitions. The 
resulting square surface is specified to have 
the conventional periodic boundary condi- 
tions with respect to 0 2 adsorption and sur- 
face reaction. Finally, the reaction is run on 
the resulting surface and a time-averaged 
steady-state rate is determined. Results for 
two of the cases we have studied are re- 
ported here: in one case diagonal actions 
were allowed, and in the other case diagonal 
actions were not allowed. In both cases 0 2 
adsorption can occur an any site pair and 
the base set of parameter values was used. 

The smallest possible base array is a 1 x 
1 array. The only configurations possible are 
the all-a surface and the all-/} surface. These 
surfaces are also possible surfaces for all 
larger base arrays and will not be mentioned 
below. Both of these surfaces are inactive 
for the two cases considered here. 

The next larger base array is a 2 x 2 array. 
Simulations were performed with surfaces 
in which the base arrays were propagated 
periodically eight times in both directions, 
with the entire 16 x 16 arrays having con- 
ventional periodic boundary conditions for 
reaction and 0 2 adsorption. All surfaces 
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Fraction of p sites - 0.3 
Rate-per-c-srte for most 
active a site = 0. 1 9 
Ave. rate-per-a-site « 0.082 


n* 





Fraction of (5 sites * 0.4 
Rate-per-a-site for most 
active a site - 0.22 
Ave. rate-per-a-site = 0.12 


Fraction of p sites « 0.5 
Rate-per-a-site for most 



Fraction of p sites = 0.6 
Rate-per-c-site for most 
active a site = 0.21 
Ave. rate-per-a-site = o.i i 


Fraction of p sites = 07 
Rate-per-a-site for most 
active a site * 0.20 
Ave. rate-per-a-site = 0.028 


Fig. 13. Sections showing the site configurations surrounding the most active a sites present in 
individual random ct-/3 surfaces for the case in which 0 : cun adsorb on any site pair and for the base 
set of parameters, a sites are represented as large black circles and /3 sites are represented as small 
black dots, with adsorbed species not shown. The most active a site is the central site in each section 
and is marked with a white dot at its center. For this set of reaction rules and parameter values, the 
rate-per-a-site for the most active a sites in random surfaces is about 0.2 at all active fraction-/?. 


with 2x2 base arrays were inactive for the 
two cases considered here, with two excep- 
tions. One exception is a 50-50 mixture of 
a and /3 sites arranged in alternating nondi- 
agonal rows of a and /3 sites. This surface 
is active for the case considered in which 
diagonal actions are not allowed, and the 
rate-per-a-site is 0.06. The rate-per-a-site 
for a 50-50 mixture of a and ft sites arranged 
randomly is 0.13. 

The other exception is a 50-50 mixture of 
a and (3 sites arranged in alternating diagonal 
rows of a and /3 sites, a configuration which 
looks like a “checkerboard” of a and (3 
sites. This surface is active for the case con- 
sidered in which diagonal actions are al- 
lowed, and the rate-per-a-site is 0.17, higher 
than for the 50-50 random surface. 

The c(2 x 2)-Sn/Pd(100) surface alloys 
studied experimentally by Logan and Paffett 
(42) have a checkboard configuration of Pt 


and Sn atoms, as determined by low energy 
electron diffraction and surface composition 
measurements. They found that the surface 
was active for CO oxidation when exposed 
to 16-TorrCO and 8-Torr0 2 at 170°C, how- 
ever, the checkerboard structure was dis- 
rupted as Sn became oxidized and presum- 
ably formed SnO x patches on top of the 
surface. An interesting question is whether 
the checkerboard Pt-Sn surface alloy would 
be stable and have a high activity for CO 
oxidation under milder conditions than used 
in (42). 

The next larger base array is a 3 x 3 array. 
This base array is sufficiently large that a 
large number of interesting surface configu- 
rations are possible. For each specified ratio 
of a and (3 sites, each distinct configuration 
of sites is a statistical “combination.” A 
search of possible statistical combinations 
of a and ft sites was performed in the follow- 
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ing manner in order to achieve a confidence 
limit of >99.9% of trying all nonequivalent 
combinations: A configuration in the base 
array which does not have two- or fourfold 
rotational symmetry is equivalent, with re- 
spect to reaction, to the three other distinct 
statistical combinations that are formed by 
rotation of this nonrotationaliy symmetric 
configuration. By running a number of trials 
with random generation of site configuration 
equal to twice the number of possible statis- 
tical combinations, at least one of these four 
equivalent configurations will be tried 
within a confidence limit of >99.9%. During 
the same number of trials, at least one of 
the four equivalent configurations of all 
combinations which do not have two- or 
fourfold rotational symmetry will also be 
tried within this confidence limit. Configu- 
rations in the square base arrays that have 
two- or fourfold rotational symmetry were 
run deliberately by specifying the site con- 
figurations manually in separate runs. 

The search for optimal site configurations 
with 3 x 3 base arrays was performed with 
surfaces in which the base arrays were prop- 
agated periodically seven times in both di- 
rections, with the entire 21 x 21 arrays 
having conventional periodic boundary 
conditions for reaction and 0 2 adsorption. 

Figures 14 and 15 summarize the results 
of the searches performed with 3x3 base 
arrays. Figure 14 is for runs in which diago- 
nal actions are not allowed and Fig. 1 5 is for 
runs in which diagonal actions are allowed. 
The notation below each pattern gives the 
fraction-/}, the rate-per-a-site, and ratc-per- 
a-site over a random surface with the same 
fraction-/}. The effect of disallowing or 
allowing diagonal actions on the most active 
site configurations is clear: in Fig. 14 the a 
and (3 sites tend to be arranged in horizontal 
and vertical rows, whereas diagonal group- 
ings of a sites and /}-/} site pairs are preva- 
lent in Fig. 15. Note the similarity between 
the local configurations of the most active 
a sites in random surfaces, shown in Fig. 
13, and the optimal 3x3 base-array config- 
urations found for the same rules and param- 


eters in Fig. 15, especially the lower left 
configuration. 

For the case in which diagonal actions are 
not allowed, all of the patterns in Fig. 14 are 
more active than the checkerboard pattern, 
w hich is inactive. For the case in which diag- 
onal actions are allowed, all of the patterns 
in Fig. 15 are less active than the checker- 
board pattern. 

In Fig. 15, one can form the optimal three- 
fi-siic pattern simply by adding one j3 site 
to the base array of the optimal two-/}-site 
pattern. Adding one /} site to the base array 
of the optimal three-/}-site pattern forms the 
optimal four-/3-site pattern, and so forth to 
form the optimal five-/}-site pattern. A simi- 
lar progression can be done in Fig. 14, ex- 
cept that movement of a (3 site is required 
to go from the optimal three-/}-site pattern 
to the optimal four-/}-site pattern. 

The results of Fig. 3 for random surfaces 
over which 0 2 can adsorb on any site pair 
can be compared with the results of Fig. 15 
for optimal 3x3 base-array patterns with 
the same 0 2 adsorption rules. In most cases 
in which the 3x3 patterns are active, they 
are somewhat more active than a random 
surface at the same fraction-/}. The 3x3 
patterns are inactive below fraction-/} = 
0.22 and above fraction-/} = 0.56, whereas 
random surfaces are active except at 
fraction-/} = 0 and at fraction-/} = 0.8 and 
above. 

The rate-per-a-site of the most active a 
site was about 0.2 reaction events per CO 
collision for all random surfaces corre- 
sponding to the diagonal-pair rules and 
fraction-/} shown in Figs. 14 and 15. This 
rate is higher than the rates found for the 
checkerboard pattern and the patterns in the 
3x3 search. This suggests that a base array 
larger than the 3x3 base array may be 
required in order to construct the overall 
optimal pattern of a and /} sites. 

The next larger base array is the 4x4 
base array. We did not do a complete study 
of 4 x 4 base arrays since the number of 
possible configurations is very large. From 
observing the patterns found in the 3x3 
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2ft (0.22) P sites 
Rat e-per-a- Site = 0.064 
Ave. rale-per-a-site for 
random surface - 0.043 




4/9 (0.44) p sites 
Rate-per-a-site * 0.1 5 
Ave. rate-per-a-site for 
random surface = 0.081 


3/9 (0.33) p sites 
Rate-per-a-site * 0.1 1 
Ave. rate-per-a-site for 
random surface - 0.069 



5/9 (0.56) P sites 
Rate-per-a-srie = 0.14 
Ave. rate-per-a-site for 
random surface = 0.067 


6/9 (0.67) p sites 
Rate-per-a-site = 0.10 
Ave. rate-per-a-sile for 
random surface = 0.032 


Fig. 14. The most active configurations found for 3x3 base arrays for cases in which diagonal 
actions are not allowed. 0 : can adsorb on any nondiagonal site pair, and the base set of parameters 
was used with the exception of not allowing diagonal actions, a sites are represented as large black 
circles and /3 sites are represented as small black dots. The case in which there is one site in the 
nine-site base array is inactive due to CO poisoning, and the cases in which there arc seven or eight 
ft sites in the base array are inactive due to oxygen poisoning of the a sites. 


base array search, we inferred that the 
4x4 base-array configuration shown in the 
upper right quadrant of Fig. 16 would have 
high activity for the case in which diagonal 
actions are allowed. The rate-per-a-site 
over this "zig-zag" pattern is higher than 
any found in the 3 x 3 base array search 
but is 4% less than over the checkerboard 
pattern shown in the upper left quadrant of 
Fig. 16. The zig-zag pattern is probably less 
active than the checkerboard because an ox- 
ygen atom adsorbed on an a site can be 
removed by CO molecules adsorbed on 
other a sites along only one diagonal direc- 
tion (i.e. , the a chains are isolated from each 
other), and a buildup of adsorbed oxygen 
atoms in one section of an a chain can lead 
to deactivation of this section of the chain. 
In the checkerboard, an oxygen atom ad- 


sorbed on an a site can be removed by CO 
molecules adsorbed on other a sites along 
two diagonal directions. 

Adding one a site to the 4 x 4 base array 
that forms the zig-zag pattern produces the 
"zig-zag + la" pattern shown in the lower 
right quadrant of the figure. This pattern 
has a lower rate-per-a-site than the pattern 
above it but has the highest overall rate of 
any surface we have identified. Note that 
the zig-zag chains of a sites in the zig- 
zag + la pattern are connected to each other 
at points every four a's along the chains, 
thus reducing the chance that a section of 
a chain will be deactivated by adsorbed oxy- 
gen. The "checkerboard + la" pattern, 
formed by' adding one a site to the 4x4 
base array that forms the checkerboard pat- 
tern, is 0.2% less active than the zig-zag + 
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2/9 (0.22) p sites 
Rate-per-a-site « 0.044 
Ave. rate-per-a-site for 
random surface « 0.050 


3/9 (0.33) p sites 
Rate-per-a-SJte - 0.12 
Ave. rate-per-a-site for 
random surface - 0.098 
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YYY 

YYY 
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4/9 (0.44) P sites 
Rate-per-a-site = 0.16 
Ave. rate-per-a-site for 
random surface = 0.13 


5/9 (0.56) p sites 
Rate-per-a-srte = 0.13 
Ave. rate-per-a-site for 
random surface = 0.12 


Fig. 15. The most active configurations found for 
3x3 base arrays for cases in which diagonal actions 
are allowed. 0 : can adsorb on any site pair, and the 
base set of parameters was used, a sites are represented 
as large black circles and /? sites arc represented as 
small black dots. The case in which there is one /3 
site in the nine-site base array is inactive due to CO 
poisoning, and the cases in w hich there are six, seven, 
or eight /? sites in the base array arc inactive due to 
oxygen poisoning of the a sites. 


la pattern. Subtracting one a site from the 
4x4 base arrays of the two top patterns in 
the figure results in decreases in both the 
rates-per-a-site and the overall rates. 

SUMMARY 

This study demonstrates how high activ- 
ity for CO oxidation can be obtained over 
a composite material composed of a highly 
interspersed mixture of one type of site, a, 
that adsorbs CO and 0 2 and another type 
of site, /?, that adsorbs 0 ; without significant 
CO inhibition. As long as 0 2 can adsorb on 
a pair of /3 sites, this mode of 0 2 adsorption 
will predominate at low temperatures, 
where the CO desorption probability is low, 
over other possible modes of 0 2 adsorption. 
The main effect of 0 2 adsorption on a-a 


or a-p pairs is partial deactivation of the 
surface at high fraction-/?, and this deactiva- 
tion will be minimal when the 0 2 sticking 
probability on a sites is low relative to the 
CO sticking probability, except at extremely 
high fraction-/?. For some reaction rules and 
parameter values, specific a-/? ratios or spe- 
cific site configurations are inactive. For 
most reaction rules and parameter values, 
a roughly 50-50 mixture of the two types of 
sites will give the greatest activity per unit 
total surface area. This result is determined 


::::: 

checkerboard 



zig-zag 



checkerboard + la zig-zag + la 


Fig. 16. 4 x 4 base-array patterns with high activity 
for the case in which 0 2 can adsorb on any site pair, 
diagonal actions are allowed, and the base set of param- 
eters is used. The two top patterns are composed of 
8/16 (0.5) /3 sites. The checkerboard pattern (also a 2 x 
2 base-array pattern) has the highest rate-per-a-silc, 
0.174, of any pattern studied in this work. The zig-zag 
pattern has a ratc-pcr-a-site of 0.167. Random surfaces 
with this fraction of /? sites have an average ratc-per- 
a-site of 0.13. The overall rates for the checkerboard, 
zig-zag, and random patterns arc 0.087, 0.084, and 
0.065, respectively. The two bottom patterns are com- 
posed of 7/16 (0.44) /? sites. The two bottom patterns 
have lower rates-per-a-sile but higher overall rates than 
the two top patterns. The rates-per-a-site for the check- 
erboard + la, the zig-zag + la, and random patterns 
with the same fraction-/? are 0.1575, 0.1580, and 0.13, 
respectively. The overall rates for the checkerboard + 
la, the zig-zag + la, and random patterns with the 
same fraction-/? are 0.08857, 0.08875, and 0.073, re- 
spectively. The zig-zag + la pattern has the highest 
overall rate of any pattern studied in this work. 
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by the reaction stoichiometry and the fact 
that the two reactants primarily adsorb sep- 
arately on the two different types of sites. 

The sites in a randomly distributed mix- 
ture of the two types of sites have widely 
differing activities which depend on the 
local site configurations. The local site 
configurations of the most active sites in a 
random surface are similar to site configu- 
rations found in the search for optimal con- 
figurations. The site configurations found in 
the search for optimal configurations were 
only about 20% more active than the random 
surfaces of the same overall composition. 
This small difference may be due to the rela- 
tively simple steric requirements of CO and 
0 2 adsorption. We expect that similar 
searches for optimal site configurations will 
be more fruitful for reactions with more 
complex steric requirements. 

More detailed experimental studies of re- 
action kinetics over composite catalysts are 
required in order to advance our under- 
standing of low-temperature CO oxidation. 
The present work emphasizes the need for 
measurements of surface composition and 
surface structure. The use of scanning nano- 
probe techniques with model catalyst sys- 
tems, such as Pd-Sn surface alloys (42), 
should be especially useful. One particularly 
interesting question is what the relative con- 
tributions are of (a) enhanced reaction at 
the perimeter of noble-metal particles, pre- 
dicted by the simulations of Brosilow ct al. 
(9) and (b) reaction over highly interspersed 
mixtures of the two components, the pres- 
ence of which have been identified experi- 
mentally by Hoflund and co-workers (39 t 
40) and the kinetics of which have been sim- 
ulated in this work. 
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Appendix I, Program Patchwork Listing, 13 pages total 


APPENDIX I 

PROGRAM: Patchwork (Microsoft QuickBasic) 

This is the program used for random surfaces. 


CLS 

CLEAR 

PRINT "Written by A jay Badlani and Dr. Richard K. Herz, Chemical" 

PRINT "Engineering, Mail Code 0310," 

PRINT "University of California at San Diego, La Jolla, CA 92093-0310," 
PRINT "Phone: (619) 534-6540, internet: rherz0ucsd.edu, bitnet:" 

PRINT "rherz0ucsd. bitnet 
PRINT "All rights reserved." 

PRINT "" 

PRINT "Hit any key to continue:" 

INPUT MX 

REM define variables 
DEFINT A-Z 
REM dimension arrays 

DIM black%{4) 'pen specification for black ovals (CO) and dots (B sites) 

DIM shade%{4) 'pen specification for shaded ovals (O atoms) 

DIM rectangle% (4 ) 'coordinates of oval in subgraphics subroutine 
DIM site% (100, 100) 'site array, see all set-up subroutines for max. 

'Ispan% and Jspan% values 

DIM rxnarray ! (100, 100) 'array matching site% that keeps track of no. of 

'reaction events on each site 
DIM SM%(10) 'for subcross subroutine 

DIM ST%(10) 'for subcross subroutine 

DIM Nrxn! (1000) 'number of reaction events in M (100 MC) time steps 
DIM VR! (100) 'used in subranO 
REM Specify Parameter Values 

yco ! =2/3 ! yco is the mole fraction of CO in the CO-02 mixture 
Mco!= 28 'Molecular weight of CO 
Mo 2!=32 'Molecular weight of 02 

y ! =1 / ( 1+ ( ( 1-yco ! ) /yco ! ) *SQR (Mco ! /Mo2 ! ) ) 'Ratio of CO collisions to the 

'total number of collisions 

REM reseed random number generator "RND" 

REM use "TIMER, " but first convert it to an integer (so can later repeat 
REM an old run exactly using same seed) 

TSEED ! = (TIMER-4 3201 !)* (655361/86401!) *65536 max integer range, 86401 

'max sec from midnight 

ISEED%=FIX (7SIE0! ) ' ISEED% should be from -32768 to +32766 

RANDOMIZE ISEED% 

REM Set up graphics by defining bit patterns for the drawing "pen" 

black% ( 0 ) =&HFFFF :black% ( 1 ) =SHFFFF : black% (2) =&HFFFF : black % ( 3 ) =&HFFFF 
shade% (0) =&H4422 : shade% ( 1 ) =&H8811 : shade! (2 ) =&H4422 : shade% (3)-&H8811 
labelOO: 'label to re-run program 

CLS 

PRINT "This is a Monte Carlo simulation of CO oxidation on a surface with" 
PRINT "two types of sites, A & B." 

PRINT "Empty B sites are marked with dots. Empty A sites are blank." 
PRINT "Black circles are CO which adsorb on a single A site only." 

PRINT "Grey circles are O atoms. 02 needs two empty sites to adsorb." 
PRINT "The input options selected determine whether actions involving" 
PRINT "next nearest neighbor sites (on diagonals) can occur and whether" 
PRINT "02 can adsorb on AB pairs, AB & BB pairs, AB & AA pairs, any pair," 
PRINT "or just B3 pairs." 
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PRINT "" 

PRINT "Rate constants are read from the file 1 Patchwork Input Data* and" 
PRINT "are equal" 

PRINT "to the probability that an event will happen when the" 

PRINT "corresponding" 

PRINT "site-occupancy configuration is chosen randomly." 

PRINT ,,n 

PRINT "Output data are written to the file * Patchwork Output ###’ where" 
PRINT "### is" 

PRINT "an index number assigned automatically. The index number is stored" 
PRINT "in the file" 

PRINT "' Patchwork Index 1 ." 
label RE AD: 


PRINT "" 

PRINT "TURNING OFF COLOR AND SCREEN SAVER SUGGESTED TO SPEED PROGRAM!!!!" 
PRINT "" 

PRINT "EDIT AND CLOSE file ’Patchwork Input Data* before continuing!" 
PRINT "" 

PRINT "Hit any key when ready to read input data file:" 

INPUT MX 
CLS 

OPEN "Patchwork Input Data" FOR INPUT AS #1 
INPUT # 1 , M ’*100 MC steps 

INPUT #l f IX ' 0=checker, l=random, 2=4Ax2B stripe, 3=cross, 4= 
’Donut, 5= read "Patch OLD sites" file 
INPUT #l,F3sire! 'Nominal fraction of B sites for IX = 1, random 

1 distribution 


INPUT # 1 , NN '0 means diagonal actions not allowed 

INPUT #1,B3 '0 for 02 ads on AB pairs, 1-AB & BB, 2=AB & AA, 3=any 

'pair, 4=3B pairs 

INPUT # 1 , kccadsorb ! 'Probability CO will adsorb on collision with 

’an empty A site (use value from 0 to 1) 

INPUT #1, kcodesorb! ’Probability that an adsorbed CO on an A site 

'will desorb when selected (use value from 0 tol) 


INPUT # 1 , kcxaisorb ! ’Probability that an 02 will adsorb on 
’collision with an empty site pair determined by BB (use value from 0 to 1) 
INPUT #l,kreact! ’trys at reaction per CO collision with an A site 

CLOSE #1 

PRINT "M (# 100 MC time steps) = ",M 

PRINT "IX (0=checker, l=random, 2=4Ax2B stripe, 3=cross, 4=Donut, 5=old sites)" 
PRINT "= ",IX 

PRINT "FBsite! (fraction of 3 sites for IX = 1, random distribution)" 
PRINT "=", FBsite! 

PRINT "NN (0 means diagonal actions not allowed) = " , NN 

PRINT "BB (0 for C2 ads on AB pairs, 1=AB & BB, 2=AB & AA, 3=any pair," 
PRINT " 4 =BB pairs) = " , BB 
PRINT "kcoadsorb! =", kcoadsorb ! 

PRINT "kcodesorb! =", kcodesorb ! 

PRINT "koxadsorb! =", koxadsorb ! 

PRINT "kreact! =",kreact! 

PRINT "" 

PRINT "For uniform, all A surface, BB needs to be = 2" 

PRINT "For IX=5, file ’Patch OLD sites' for 1=42, J=26 array must be" 
PRINT " present 
PRINT "" 

PRINT "Hit 1 to re-read input data, any other key to continue:" 

INPUT MX 

IF MX=1 THEN GOTO labelREAD 
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REM **** RE-activate the next 7 lines by removing leading {') to do SCAN ***** 
' IX=1 'override input params to get random distribution 
' BB=2 'override input params on BB, CHECK TO MAKE SURE THIS IS WHAT 
' YOU WANT !!!!!! 

' M=1 5 'override input params on M, CHECK TO MAKE SURE THIS IS WHAT 
1 YOU WANT !!!!!!! 

'FOR ZII=1 TO 9 ***** TEST (should be 9) **** 

' FBsite ! = ( . 1*ZII) ’scan range of FBsite! 

REM **** repeat 10 times FOR each value of FBsite! *** 

'FOR ZI-1 TO 10 '**** TEST (should be 10) **** 

REM *********************************************************** 

REM read index file for output file names 
OPEN "Patchwork Index" FOR INPUT AS #1 
INPUT #1, indexA 
CLOSE #1 

REM update index file 

OPEN "Patchwork Index" FOR OUTPUT AS #1 

indexA=indexA+l 

PRINT #1, indexA 

CLOSE #1 

REM open output file 

fileout$="Patchwork Output "+STR$ (indexA) 

OPEN f ileout$ FOR OUTPUT AS #1 

PRINT #1, "Patchwork Output ", STR$ (indexA) 

PRINT #1,"" 

PRINT #1, "Initial state (0=checker, l=random, 2=4Ax2B stripe, 3=cross, 
4=Donut, 5=old sites) = ",IX 
PRINT #1,"" 

PRINT #1, "FBsite! (fraction of B sites for IX = 1, random distribution)" 
PRINT #1,"= ", FBsite! 

PRINT #1,"" 

PRINT #1,"NN (0 means diagonal actions not allowed) = ",NN 
PRINT #1,"" 

PRINT #1 , "BB (0 for 02 ads on AB pairs, 1=AB & BB, 2=AB & AA, 3=any pair," 
PRINT # 1 , " 4=BB pairs) = ",B3 
PRINT #1,"" 

PRINT #1, "kcoadsorb! = ",kcoadsorb! 

PRINT #1, "kcodescrb! = ",kcodesorb! 

PRINT #1, "koxadsorb! = ",koxadsorb! 

PRINT #l,"kreact! = ",kreact! 

PRINT #1, "yco!=",yco! 

PRINT #1,"" 

PRINT #1 , " I SEED% =" , ISEED% 

PRINT #1,"" 

REM set-up initial state 

IF IX=0 THEN G0SU3 subchecker 
IF IX=1 THEN G0SU3 subrandom 
IF IX=2 THEN G0SU3 substripes 
IF IX=3 THEN G0SU3 subcross 
IF IX=4 THEN G0SU3 subdonut 
IF IX=5 THEN G0SU3 subsitefile 
REM return from initial state set-up subroutines 

REM calculate NStotal! the following way to not exceed integer limit 
NStotal ! =Ispanl 
NStotal ! =NStotal ! * Jspan% 

FAactual ! =NAsite ! /NStotal ! 

FBactual ! =l-FAactual ! 

REM krxn! is the number of trys at reaction to the total trys 


24 



Appendix I, Program Patchwork Listing, 13 pages total 


krxn ! = (y ! *kreact ! ) / (1+ (y ! *kreact ! ) ) 

PRINT #1, "The total number of sites =",NStotal! 

PRINT #l,"The number of A sites = ,, ,NAsite! 

PRINT #1, "FAactual! = ",FAactual! 

PRINT #1, "FBactual ! (compare to FBsite! for IX = 1) = ",FBactual! 

PRINT #1 , "krxn ! =",krxn! 

PRINT #1,"" 

PRINT #l,"Rate is no. of reaction events per CO collision with an A site, 
PRINT #1, "averaged over 100*NStotal! trys (times through loop)." 

PRINT #l, H Dtime! is the no. of CO collisions per A site from time - 0." 
PRINT #1,"" 

PRINT #1, "Dtime! ", ", ", "Rate" 

REM display initial state 
CLS 

GOSUB subgraphics 

REM initialize reaction event counter array 
FOR 1=1 TO Ispan% 

FOR J=1 TO Jspan% 

rxnarray ! ( I , J ) — 0 ! 

NEXT J 
NEXT I 

REM start main iteration loop 
ON BREAK GOSUB labelBREAK : BREAK ON 
FOR KK=1 TO M 

REM initialize reaction event counter 
Nrxn ! (KK) =0 

REM Set flag IFF so "RND" will be re-seeded in subranO every 10 KK (every 
REM 1000 MC steps) on average 

REM in order to increase repetition period of "RND" 

GOSUB subranO 

IF RANO ! < . 1 THEN IFF=0 

FOR K=1 TO 100 '****** TEST (should be 100) ********* 

FOR LA=1 TO Ispan% 'split into two loops so don't exceed integer limit with 

'large Ispan% * Jspan% 

FOR LB=1 TO Jspan% 

REM Check for a reaction krxn! fraction of the time and for an adsorpticr.- 
REM desorption event (1-krxn!) of the time. 

GOSUB subranO 

IF RAN0!>krxn! THEN GOTO label88 'then look at an adsorption event 
REM check to see if reaction happens 
GOSUB subpairpick 

REM REACT if 0 on A site and CO on A site 

IF site% (x, y ) =2 AND site% (xn% , yn%) =1 THEN 
GOSUB subrxn 
GOTO labellO 
END IF 

REM REACT if CO on A site and O on A site 

IF site% (x, y ) =1 AND site! (xn% , yn% ) =2 THEN 
GOSUB subrxn 
GOTO labellO 
END IF 

REM REACT if CO on A site and O on B site 
IF site% (x,y) =1 AND site% (xn%, yn% ) =12 THEN 
G0SU3 subrxn 
GOTO labellO 
END IF 

REM REACT if O on B site and CO on A site 

IF site% (x, y) =12 AND site% (xn% , yn% ) =1 THEN 
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GOSUB subrxn 
GOTO label 10 
END IF 

GOTO labellO 
label88 : 

REM look at a CO event y! fraction of the time and at an oxygen event 

REM (1-y!) of the time 

GOSUB subranO 

IF RANO ! >y ! THEN GOTO label888 'then go look at an oxygen event 

REM check for CO adsorption or desorption 
REM pick random site, site%(x,y) 

GOSUB subranO 

x=INT( ( Ispan%*RANO ! ) +1 ! ) 

GOSUB subranO 

y=INT ( ( Jspan%*RANO ! ) +1 ! ) 

IF site% (x, y ) =0 THEN 
G0SU3 subcoadsorb 
GOTO labellO 
END IF 

IF site% (x, y ) =1 THEN GOSUB subcodesorb 

GOTO labellO 
label888 : 

REM check to see of 02 adsorbs 

GOSUB subpairpick 

REM ADSORB 0 if B3<4 and empty A site and empty B site 
IF BB<4 AND site%<x,y)=0 AND site% (xn%, yn% ) =1 0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if B3<4 and empty B site and empty A site 
IF BB<4 AND site% (x, y ) =1 0 AND site% (xn% , yn% ) =0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if BB=1 and empty B site and empty B site 
IF BB=1 AND site% (x,y)-10 AND site% (xn% , yn% ) -10 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if 33=2 and empty A site and empty A site 
IF BB=2 AND site! (x, y ) =0 AND site% (xn% , yn% ) =0 THEN 
G0SU3 suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if B3=3 and empty A site and empty A site 
IF BB=3 AND site%(x,y)=0 AND site% (xn% , yn% ) =0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if B3=3 and empty B site and empty B site 
IF BB=3 AND site% (x , y ) =1 0 AND site% (xn%, yn%) =10 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if EB=4 and empty B site and empty B site 

IF BB=4 AND site% (x, y ) =1 0 AND site% (xn% , yn% ) =1 0 THEN 
G0SU3 suboxadsorb 
GOTO labellO 
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END IF 

labellO: 

'GOSUB subgraphics **** activate by removing leading {') to check logic 

’changes by displaying each trial ****** 

NEXT LB 
NEXT LA 
NEXT K 

GOSUB subgraphics 

IF NAsite ! =0 THEN GOTO labelSKIP 

REM at this point Nrxn ! (KK) equals no. of reaction events in 1 00 ! *NStotal ! 
REM total times through loop. 

Nrxn! (KK) =Nrxn ! (KK) / (100 ! *NStotal ! * (1-krxn ! ) *y ! *FAactual ! ) 

REM Nrxn! (KK) now equals no. of reaction events per CO collision with an A 
REM site 

Dtime ! =KK*100* (1-krxn !) *y ! 'Dtime! is total no. of CO collisions per A 

'site since time - 0 

labelSKIP: 

LOCATE 21 

PRINT USING "### f .#li AAAAI, ;KK,Nrxn! (KK) 

PRINT #1, USING # . ### ;Dtime ! , Nrxn ! (KK) 

NEXT KK 
labelBREAK : 

REM If allowed to complete FOR-NEXT loop, KK becomes M+l, so reset it to 
REM equal M. 

REM If BREAK sent, don't use incomplete series in average below so reset 
REM anyway. 

KK=KK-1 

REM **** DE-activate the next 15 lines by adding leading (') to do SCAN **** 
LOCATE 22 

'PRINT "hit any key to re-display final configuration:" 

'INPUT MX 

GOSUB subgraphics 
LOCATE 21 

PRINT USING "#*# #.### A ; KK, Nrxn ! (M) 

PRINT " 

'PRINT "hit ar.y key to continue: " 

'INPUT MX 

IF KK>1 5 THEN MS=KK-14 ELSE MS=1 
LOCATE 1 

PRINT "Rate 0 M (# 1 00 *NStotal ! trys) " , MS, "to" , KK 
FOR KS=MS TO KK 

PRINT USING "# . ### AAAA ”;Nrxn! (KS) 

NEXT KS 

★★*****VtT*:**'******icx>t******x * *★★*****★****★*★*•*★*★***★*★*'*★ 

REM ave the last half of the run to get an average rate 

Avestart %=INT ( . 5*KK) 'start averaging 50% through series - can change 

'factor if desired 

Nrxn ! (M+l ) =0 'place holder for the average of (KK+l-Avestart ) Nrxn! to 
'average together 
FOR KS=Avestart % TO KK 

Nrxn! (M+l ) =Nrxn! (M+l ) +Nrxn ! (KS) / (KK+1-Avestart% ) 

NEXT KS 
LOCATE 22 

PRINT "Rate =" , Nrxn ! (M+l ) 

GOSUB subanalysis 

REM finish writing to output file 

PRINT #l,"The following is a ' snapshot * of the final site configuration," 
PRINT #1, "site% (I, J) , " 
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PRINT #1," listed as all Jspan% J's for 1=1, then all Jspan% J's for 1=2," 
PRINT #1, "etc. , up to I=Ispan%" 

PRINT #1, "Ispan% = ", Ispan% 

PRINT #1, "Jspan% =", Jspan% 

FOR 1=1 TO Ispan% 

FOR J=1 TO Jspan% 

PRINT #1, site% (I, J) 

NEXT J 
NEXT I 
PRINT #1,"" 

PRINT #l,"The following is the reaction event counter array" 

PRINT #1, "rxnarray ! (I, J) , " 

PRINT #1," listed as all Jspan% J's for 1=1, then all Jspan% J's for 1=2," 
PRINT #1, "etc. , up to I=Ispan% : " 

FOR 1=1 TO Ispan% 

FOR J=1 TO Jspanl 

PRINT #1, rxnarray ! (I, J) 

NEXT J 
NEXT I 
PRINT #1,"" 

PRINT #1 , "END OF DATA" 

CLOSE #1 

REM update summary file 

OPEN "Patchwork Summary" FOR APPEND AS #1 

PRINT #1, USING "*#*_,######_,_ #_,_ #.###_,_ #_,_ #.### 

#.###_,_ #.###_,_ #.###_,_ 

#.###"; indexA, ISEEDI, IX, FBactual ! , NN, BB, Nrxn ! (M+l ) , kcoadsorb! , kcodesorb ! , kcxad 
sorb ! , kreact ! 

CLOSE #1 

REM **** RE-activate r.ext 2 lir.es by removing leading (') to do a SCAN 
***************** 

1 NEXT ZI 
1 NEXT ZII 

REM **** DE-activate next 7 lines by adding leading (') to do a SCAN 
******************* 

LOCATE 21 
PRINT " 

PRINT 11 

PRINT "ENTER 1 TO RERUN: 

PRINT " 

INPUT MX 

IF MX=1 GOTO labelOO 

REM ***************** xx ******** * ************************** 

BREAK OFF 
END 

subranO : 

REM This routine is similar "RANO" in the book "Numerical Recipes (FORTRAN) by 
REM W. H. Press, et al." 

REM It "scrambles" the "RND" numbers to cut down on pair correlations. 

REM I have also added a random re-seeding of "RND" to increase its repetition 
REM period 

REM from originally on the order of 5E07 calls (11,000 MC steps for 1092 
REM array) . 

IF IFF=0 THEN 
IFF—l 

REM re-seed "RND" randomly with random seed to to increase its 
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REM repetition period 
ZR% =32767* END 
IF RNDC.5 THEN ZR%=-ZR% 

RANDOMIZE ZR% 

REM the next lines "set up" the scrambling array for RANO 
FOR JR-1 TO 97 
dum! =RND 
NEXT JR 
FOR JR=1 TO 97 
VR! (JR)=RND 
NEXT JR 
YR! =RND 

END IF 

JR=1+INT (97 ! * YR ! ) 

YR!=VR! (JR) 

RANO !=YR! 

VR! (JR) =RND 

RETURN 
subchecker : 

REM initialize for A-3 checker board initial state 

Ispan%=42 ? can be changed for different array size but needs to be 
* even number 

Jspan%=26 'can be changed for different array size but needs to be 
'even number 
NAsite ! =Ispan*/2 

NAsite ! =NAsire ! x Jspan% 'calculate this way so don't exceed integer 

'limit during calculation 
FOR 1=1 TO 21 '1 to (Ispan%/2) 

10 = ( 2 * 1) -1 
IE= (2*1) 

FOR J=1 TO 13 '1 to { Jspan%/2 ) 

J0= (2 x J) -1 
JE= (2 X J) 

site* (10, JO) =0 'site is empty A site 
site* (IE, JE) =0 ’site is empty A site 
site* (10, JE) =1 0 'site is empty B site 
site* (IE, JO) =10 'site is empty B site 
NEXT J 
NEXT I 

RETURN 

subrandom: 

REM initialize for empty random A-B site distribution 
NAsite !=0! 

Ispan%=42 'can be changed for different array size 
Jspan%=26 r can be changed for different array size 
FOR 1=1 TO Ispanl 

FOR J=1 TO Jspan* 

G0SU3 subranO 

IF RANG ! >FBsite ! THEN 

site%(I,J)=0 'site is empty A site 

NAsite ! =NAsite ! +1 ! 'number of A sites to use in 

'calculating rate per A site 

GOTO label SRO 
END IE 

site! (I, J) =10 'site is empty B site 

labelSRO: 

NEXT J 
NEXT I 
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RETURN 
substripes : 

REM initialize for stripes of 4 A sites wide followed by 2 B sites wide 
Ispan%=42 ’can be changed for different array size but must be 
’multiple of repeat period of 6 
Jspan%=26 'can be changed for different array size 
NAsite ! =Jspanl* (4*Ispan%/6) 

FOR 1=0 TO 6 ’0 to ( (Ispan%/6) -1) 

FOR J=1 TO Jspan% 

site% { (1*6) +1, J) =0 'site is empty A site 

site* ( (1*6) +2, J) =0 'site is empty A site 

site% ( (I*6)+3,J)=0 'site is empty A site 

site* ( (1*6) +4, J) =0 'site is empty A site 

site% ( (I*6)+5,J)=10 'site is empty B site 

site* ( (1*6) +6, J) =10 'site is empty B site 

NEXT J 
NEXT I 
RETURN 
subcross : 

REM initialize for staggered crosses of A's 

Ispan%=40 'car. be changed but must be multiple of repeat period of 8 

Jspan%-24 'can be changed but must be multiple of repeat period of 8 

NAsite!=600 * 6 C 0 is for Ispan%=40, Jspan%=24 
' (= Ispanl * Jspan% * 1 0/1 6 ????) 

REM initialize arrays 
SM% (1) =0 
SM% (2) =0 
SM% (3) =0 
SM% (4) =10 
SM% (5) =0 
SM% (6) =10 
SM% (7) -3 
SM% (8) =10 
ST% (1) =1 
ST% (2) =4 
ST% (3) =7 
ST% (4) =2 
ST% (5) =5 
ST% (6) =8 
ST% (7) =3 
ST% (8) =6 

FOR J=0 TO 2 '0 to { Jspan%/8 ) -1 

FOR JX=1 TO 8 

FOR 1=0 TO 4 '0 to (Ispan%/8) -1) 

FOR IX=1 TO 8 

CX=ST% (JX) + IX 

IF CX>8 THEN CX=CX-8 

site% { (I*8+IX) , ( J*8+JX) ) =SM% (CX) 

NEXT IX 
NEXT I 
NEXT JX 
NEXT J 
RETURN 
subdonut : 

REM initialize for A-B checker board initial state 
Ispan%=42 'can be changed for different array size 
Jspan%=26 'can be changed for different array size 
NAsite 1 =Ispan%/2 


30 



Appendix I , Program Patchwork Listing, 13 pages total 


NAsite ! =NAsite ! *Jspan% ’calculate this way so that it doesnot exceed 

’integer limit during calculation 
FOR 1=1 TO 21 '1 to Ispan%/2 
FOR J=1 TO 26 

site%(I,J)=0 * site is empty A site 
NEXT J 
NEXT I 

FOR 1= 22 TO 42 ’Ispan%/2 to Ispan% 

FOR J= 1 TO 26 

site% { I , J) =10 ’site is a empty B site 
NEXT J 
NEXT I 
RETURN 
subsitef ile : 

REM This is set up to read the array ’’snapshot” of a previous run for an 
REM 1 = 42, J=2 6 array. 

REM The file "Patch OLD sites" must be present. 

Ispan%=42 'can be changed for different array size 
Jspan%=26 'can be changed for different array size 
NAsite ! =0 

OPEN "Patch OLD sites” FOR INPUT AS #2 
FOR 1=1 TO 42 
FOR J=1 TO 26 

INPUT# 2 , site% ( I , J) 

IF site% (I, J) =0 THEN NAsite !=NAsite!+l 
IF site% (I, J) =1 THEN site% (I, J) =0 : NAsite ! =NAsite ! +1 
IF site! (I, J) =2 THEN site %< I , J) =0 : NAsite ! =NAsite ! +1 
IF site% (I, J) =12 THEN site% (I, J)=10 
NEXT J 
NEXT I 
CLOSE #2 
RETURN 
subpairpick : 

REM this picks a random pair of sites 
REM pick random site, site%(x,y) 

GOSUB subranO 

x=INT { (Ispan%*RAN0 \ ) +1 ! ) 

GOSUB suoranO 

y=INT ( { Jspan%*RAN0 ! ) +1 ! ) 

REM pick neighboring site 

IF NN=0 THEN GOTO labelDD 

REM the following 7 lines are for: diagonal neighbor actions are OK 

labelNN : 

GOSUB subrar.0 
xn%=x+INT (3 *RAN0 ! ) -1 
GOSUB subranO 
yn%=y+IN? (3*RAN0 ! ) -1 

IF xn%=x AND yn%=y THEN GOTO labelNN 'neighbor can't be itself 
GOTO labelS99 

labelDD: 

REM the following 18 lines are for: diagonal actions not allowed 

GOSUB subranO 

IF RAN0K.5 THEN GOTO label99 
xn%=x 

GOSUB subranO 
IF RAN 0 ! < . 5 THEN 
yn%=y+l 

ELSE 
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yn%“y-l 
END IF 

GOTO Iabei999 

label99: 

yn!=y 

GOSUB subranO 
IF RAN 0 ! < . 5 THEN 
xn!=x+l 

ELSE 

xn!=x-l 
END IF 
label999 : 

REM Set indices correctly for neighbors over array boundary. 

REM This corresponds to the surface as a torus or as a semi-infinite 
REM flat 

REM surface formed of a periodic array of the sub-array computed here. 
IF xn! = 0 THEN xn%=Ispan% 

IF xn%>I spar.% THEN xn!=l 
IF yn!~0 THEN yn%=Jspan% 

IF yn%>Jspan% THEN yn%»l 

RETURN 
subcoadsorb : 

GOSUB subranO 

IF RANO ! <kcoadsorb ! THEN site%(x,y)=l 
RETURN 
subcodesorb : 

GOSUB subranO 

IF RANO ! <kcodesorb! THEN site%(x,y)=0 
RETURN 
suboxadsorb : 

GOSUB subranO 

IF RANO ! <koxadsorb ! THEN 

site% (x,y)=sice!(x,y)+2 
site% (xn%, yn! ) =site% (xn% , yn% ) +2 
END IF 
RETURN 
subrxn : 

GOSUB subranO 

IF RANO !>=. 4052345* THEN GOTO labelRXN 
Nrxn ! ( KK ) =Nrxr. ! ( KK ) + 1 ! 

IF site%(x,y)=l THEN site!{x,y)=0 
IF site%{x f v)=2 THEN site!(x,y)=0 
IF site% (x, y) =12 THEN site! (x , y ) =1 0 
IF KK>3 6 THEN rxnarray ! <x, y ) =l + rxnarray ! (x, y ) 

IF site! (xn%, yn% ) =1 THEN site% (xn% , yn% ) =0 
IF site% (xn! , yr.% ) =2 THEN site! (xn% , yn% ) =0 
IF site% (xn% , yr.% ) =12 THEN site% (xn% , yn! ) =1 0 
IF KK>36 THEN rxnarray! (xn! , yn% ) =l+rxnarray ! (xn%,yn%) 
labelRXN: 

RETURN 
subgraphics : 

FOR 1=1 TO I span! 

FOR J=1 TO Jspan! 

REM set coordinates for "ovals" (circles here) in PAINTOVAL calls 
rectangle! (0) * ( J-l) *12 
rectangle! (1) = (I— 1 ) *12 
rectangle! (2) =J*12 
rectangle! (3) =1*12 
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dot% (0)=((J-l)*12)+5 
dot% (l)=((I-l)*12)+5 
dot% (2) = ( J*12) -5 
dot% (3) = (1*12) -5 
REM update site graphics 

IF site%(I / J)=0 THEN 'site is empty A site 
CALL ZRASEOVAL (VARPTR ( rect angle% ( 0 ) ) ) 

GOTO labelGRl 
END IF 

IF site%(I,J)=l THEN 'site is CO ads on A site 
CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR (rectangle% (0) ) ) 

GOTO labelGRl 
END IF 

IF site%(I,J)=2 THEN 'site is 0 ads on A site 
CALL PEN? AT (VARPTR (shade% (0) ) ) 

CALL PAINTOVAL (VARPTR ( rectangle% { 0 ) ) ) 

GOTO labelGRl 
END IF 

IF site! (I, J) =10 THEN 'site is empty B site 
CALL ERASEOVAL (VARPTR (rectangle! ( 0 ) ) ) 

CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR (dot % (0) ) ) 

GOTO labelGRl 
END IF 

IF site% ( I , J) -12 THEN 'site is O ads on B site 
CALL PENPAT (VARPTR { shade% (0) ) ) 

CALL PAINTOVAL (VARPTR ( rect angle% (0) ) ) 

CALL PENPAT (VARPTR <black% (0) ) ) 

CALL PAINTOVAL (VARPTR (dot % (0) ) ) 

GOTO labelGRl 
END IF 

PRINT "error: in subgraphics with unknown site configuration" 

'shouldn't reach this point 
labelGRl : 

NEXT J 
NEXT I 
RETURN 
subanalysis : 

REM determine fractional coverages and statistics 
sumco ! =0 
sumox ! =0 

FOR 1=1 TO Ispan% 

FOR J=1 TO Jspan% 

IF site%(I,J)=l THEN sumco ! =sumco ! +1 ! 

IF site! (I , J) =2 THEN sumox!=sumox!+l! 

NEXT J 
NEXT I 

IF NAsite ! =0 THEN 
thetaco ! =359 ! 
thetaox ! =939 ! 

GOTO label ANA LI 
END IF 

thetaco ! =sumco ! /NAsite ! 
thetaox ! =sumcx ! /NAsite ! 
labelANALl : 

PRINT #1, "Thetas on A sites ONLY of 'snapshot* of final 
configuration: " 
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PRINT #1,"" 

PRINT #1, USING "_t_h_e_t_a C_0 =_ #.####"; thetaco ! 

PRINT #1, USING ”_t_h_e_t_a 0_X =_ #.####"; thetaox ! 

RETURN 
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APPENDIX II 

PROGRAM: Opt? (?) /3x3x7 (Microsoft QuickBasic) 

This is the program used for searching for optimal surface structures. 


CLS 

CLEAR 

PRINT "Written by Richard K. Herz, Chemical Engineering, Mail Code 0310," 
PRINT "University of California at San Diego, La Jolla, CA 92093-0310," 
PRINT "Phone: (619) 534-6540, internet: rherz@ucsd.edu, bitnet:" 

PRINT "rherz@ucsd. bitnet " 

PRINT "All rights reserved." 

PRINT "" 

PRINT "Hit any key to continue:" 

INPUT MX 

REM define variables 
DEFINT A-Z 
REM dimension arrays 

DIM black% (4 ) 1 pen specification for black ovals (CO) and dots (B sites) 

DIM shade% (4 ) 'per. specification for shaded ovals (O atoms) 

DIM rectangle% (4) 'coordinates of oval in subgraphics subroutine 
DIM site% (100, 100) 'site array, see all set-up subroutines for max. 

Ispan% and Jspan% values 

DIM maxsite% (100, ICO) 'site array to hold optimal configuration 
DIM siteO % (100,100) 

DIM Nrxn ! (1000) 'number of reaction events in M (100 MC) time steps 
DIM VR! (100) 'used in subranO 
REM Specify Parameter Values 

yco!=2/3 'yco is the mole fraction of CO in the CO-02 mixture 
Mco!=28 'Molecular weight cf CO 
Mo2!=32 'Molecular weight cf 02 
y ! =1 / (1+ ( (1-yco ! ) /yco! ) *SQR (Mco ! /Mo2 ! ) ) 

REM y! « Ratio of CO collisions to the total number of collisions 
REM reseed random number generator "RND" 

REM use "TIMER," but first convert it to an integer 

REM (so can later repeat an old run exactly using same seed) 

TSEED ! - (TIMER-43201 !)*(65536!/86401!) 

REM 63536 max integer range, 86401 max sec from midnight 
ISEED%=FIX (TSEED! ) 'ISEED% should be from -32768 to +32766 
RANDOMIZE ISEED% 

REM Set up graphics by defining bit patterns for the drawing "pen" 

black% ( 0 ) =&HFFFF : black% ( 1 ) = S HFFFF : black% (2 ) =&HFFFF : black % (3) -&HFFFF 
shade% (0 ) =&H4422 : shade% (1 ) =SH8 811 : shade% (2) =&H4 422 : shade % (3)=&H8 811 
CLS 

PRINT "This is a Monte Carlo simulation of CO oxidation on a surface with" 
PRINT "two types of sites, A & B." 

PRINT "Empty B sites are marked with dots. Empty A sites are blank." 
PRINT "Black circles are CG which adsorb on a single A site only." 

PRINT "Grey circles are O atoms. 02 needs two empty sites to adsorb." 
PRINT "The input options selected determine whether actions involving" 
PRINT "next nearest neighbor sites (on diagonals) can occur and whether" 
PRINT "02 can adsorb on A3 pairs, AB & BB pairs, AB & AA pairs, any pair," 
PRINT "or> just 33 pairs." 

PRINT "" 

PRINT "" 

PRINT "TURNING OFF MONITOR COLOR AND SCREEN SAVER IS SUGGESTED TO SPEED UP 
PROGRAM! ! ! ! " 
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PRINT "" 

PRINT "Hit any key when ready to read input data:" 

INPUT MX 

CLS 


labelASKOl : 

PRINT "Enter no. of A sites in base array (NAbase%):" 

INPUT NAbase% 

PRINT "Enter NN (0 diag no allowed, 1 diag OK):" 

INPUT NN 

PRINT "NAbase! =",NAbase% 

PRINT "NN = ",NN 

PRINT "Enter 2 if you want to revise, any other key to continue:" 
INPUT MX 

IF MX=2 GOTO labelASKOl 


M=100 'number of 100 MC steps performed 

avestart%=10 ’number of runs from start to skip in getting average rate 

IX=1 

BB=3 

kcoadsorb! =1 
kcodesorb ! =0 
koxadsorb ! =1 


kreact ! =1 


I0span%=3 ' size of case array in IOspan direction 
J0span%=3 ’size of base array in JOspan direction 

expander%=7 'number of times base array repeated in each direction 
REM krxn! is the number of trys at reaction to the total trys 
krxn !=(y! *kreact!) / (1+y ! * kreact ! ) 


REM set up name of output file 
Mtitle$=CHR$ (NAbase!+4 8 ) 

Nt itle$=CHR$ (NN+4 3 ) 

Ititle$=CHR$ (lOspar.%+48) 

Jt itle$=CHR$ ( J0s?ar.! + 4 8 ) 

Etitle$=CHR$ (expar.der%+48) 

f ileout$="Opt "+Mtitle$+" ("+Ktitle$+") / "+Ititle$+"x"+Jtitle$+"x"+Etit le$+” 
REM print to screen 


PRINT "M (# 100 MC time steps) = ",M 
PRINT "avestart! = ",avestart% 

PRINT "IX (1 = random.) = ",IX 

PRINT "NN (0 means diagonal actions not allowed) = ",NN 

PRINT "BB {0 for C2 ads on AB pairs, 1=AB & BB, 2=AB & AA, 3 = any pair, 

PRINT "4=BB pairs) = " , BB 

PRINT "kcoadsorb! =", kcoadsorb ! 

PRINT "kcodesorb! =", kcodesorb ! 

PRINT "koxadsorb! =", koxadsorb ! 

PRINT "kreact! kreact! 

PRINT "" 

PRINT "I0span%, J0span% = ", 1 0span%, J0span% 

PRINT "NAbase% = n ,NAbase% 

PRINT "expander! = ", expander! 

PRINT "" 


labelASKOO: 

PRINT "Enter 1 to start with new random array, 2 to read old file," 
PRINT "3 to input manually:" 

INPUT MZ N 

IF MZ<1 OR MZ>3 THEN PRINT "TRY AGAIN": GOTO labelASKOO 
REM set-up initial state 

IF MZ=1 THEN GOSU3 subrandom 
IF MZ=2 THEN G0SU3 suboldfile 
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IF MZ-3 THEN G0SU3 submanual 
GOSUB subexpand 

REM calculate NStotal! the following way to not exceed integer limit 
REM during calculation 
NStotal ! =Ispan% 

NStotal ! =NStotal ! * Jspan% 

FAactual ! =NAsite ! /NStotal ! 

FBactual ! =1-F Aactual ! 

REM initialize parameters and write to ouput file so run type recorded if all 
REM poison 

maxrate ! =0 
runnum! =0 
GOSUB suboutput 
label START: 

OPEN "OPT Index" FOR INPUT AS #1 
INPUT #1 , indexA 
CLOSE #1 

REM update index file 

OPEN "OPT Index" FOR OUTPUT AS #1 
indexA=indexA+ 1 
PRINT #1, indexA 
CLOSE #1 

runnum! =runnum! +1 
IF runnum! >1 THEN 

GOSUB subrandcn 'generate new random surface for next run 
GOSUB subexpand 

REM calculate NStotal! the following way to not exceed integer limit 
REM during calculation 
NStotal ! =1 scan% 

NStotal ! -NStotal ! *Jspan% 

FAactual ! =NAsite ! /NStotal ! 

FBactual ! =l-FAactual ! 

END IF 
CLS 

GOSUB subgraphics 
REM start main iteration loop 

ON BREAK GOSUB labelBREAK : BREAK ON 
FOR KK=1 TO M 

REM initialize reaction event counter 
Nrxn ! (KK) =0 

REM Set flag IFF so "RND" will be re-seeded in subranO every 10 KK (every 
REM 1000 MC steps) on average 

REM in order to increase repetition period of "RND" 

GOSUB subranO 
IF RAN0 ! < . 1 THEN IFF=0 
FOR K=1 TO 100 

FOR LA=1 TO Ispan% 'split so don't exceed integer limit in Ispan%* Jspan% 

FOR LB=1 TO Jspan% 

REM Check for a reaction krxn! fraction of the time and for an adsorption- 
REM desorption event (1-krxn!) of the time. 

GOSUB subranO 

IF RANQ!>krxn! THEN GOTO label88 'then look at an adsorption event 
REM check to see if reaction happens 
GOSUB 1 subpairpick 

REM REACT if O on A site and CO on A site 

IF site% (x, y) =2 AXD site% (xn%, yn%) =1 THEN 
GOSUB subrxn 
GOTO labellO 
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END IF 

REM REACT if CO on A site and O on A site 

IF site% (x,v) =1 AND site% (xn%, yn! ) =2 THEN 
GO SUB subrxn 
GOTO label 10 
END IF 

REM REACT if CO on A site and 0 on B site 
IF site%(x,y)=l AND site% (xn%, yn% ) =12 THEN 
GOSUB subrxn 
GOTO label 10 
END IF 

REM REACT if 0 on B site and CO on A site 

IF site% (x, y ) =12 AND si te% (xn% , yn% ) =1 THEN 
GOSUB subrxn 
GOTO label 10 
END IF 

GOTO labellO 
label8 8 : 

REM look at a CO event y! fraction of the time and at an oxygen event ( 

REM y!) of the tire 

GOSUB subranO 

IF RANO ! >y ! THEN GOTO label888 'then go look at an oxygen event 

REM check for CO adsorption or desorption 
REM pick random site, site%(x,y) 

GOSUB subranO 

x=INC ( ( I span %* RANO ! ) +1 ! ) 

GOSUB subranO 

y=INT ( ( Jspan%*RAN0 ! ) +1 ! ) 

IF si tel y) =0 THEN 
GOSUB subcoadsorb 
GOTO labellO 
END IF 

IF site! y) =1 THEN GOSUB subcodesorb 

GOTO labellO 
label8 88 : 

REM check to see zz 02 adsorbs 

GOSUB subpairpick 

REM ADSORB O if B3<4 and empty A site and empty B site 
IF BB<4 AND site! (x, y ) =0 AND site! (xn% , yn% ) =1 0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if 33<4 and empty B site and empty A site 
IF BB<4 AND site! (x, y ) =1 0 AND site! (xn% , yn% ) =0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if 53=1 and empty B site and empty B site 

IF BB=1 AND site! (x, y) =10 AND site! <xn% , yn% ) =1 0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if 3B=2 and empty A site and empty A site 
IF BB=2 AND site! (x, y ) =0 AND site! (xn! , yn! ) =0 THEN 
G0SU3 suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB O if 3B=3 and empty A site and empty A site 
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IF BB=3 AND site%(x,y)=0 AND site% (xn%, yn% ) -0 THEN 
GO SUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if BB=3 and empty B site and empty B site 
IF BB=3 AND site% (x, y) =10 AND site% (xn% , yn% ) =1 0 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

REM ADSORB 0 if BB=4 and empty B site and empty B site 

IF BB=4 AND si te % (x , y ) =1 0 AND site% (xn%, yn%) =10 THEN 
GOSUB suboxadsorb 
GOTO labellO 
END IF 

labellO: 

'GOSUB subgraphics * activate by removing leading (') to check logic 
NEXT LB 
NEXT LA 
NEXT K 

GOSUB subgraphics 

IF NAsite ! =0 THEN PRINT "NAsite ! « 0 !!!!!!!!!!!!!!!": GOTO labelSKIP 
REM at this point Nrxn ! (KK) equals no. of reaction events in 100 ! *NStota 
REM total times through loop. 

Nrxn ! (KK) =Nrxn ! (KK) / (100 ! * (1-krxn ! ) *y ! * NAsite ! ) 

REM Nrxn! (KK) now equals (no. of reaction events per A site) per (CO 
REM collision per A site) 

Dtime ! =KK*1 00* (1-krxn! ) *y! 

REM Dtime! is total number CO collisions per A site since time = 0 
labelSKIP: 

LOCATE 21 
PRINT f ileout $ 

PRINT USING **#£# # . ###~ ~^ A ; KK, Nrxn ! (KK) 

REM ******** TERMINATE RUNNING ON DEAD SURFACE ************** 

IF Nrxn! (KK)=0 AND kcodesorb ! =0 THEN GOTO labelSTART 
IF Nrxn! (KK) =0 THEN 
killcheck%=0 
FOR 1=1 TO Is?an% 

FOR J=1 TO Jspanl 

IF sitel(I,J)=0 OR site% (I , J) =1 THEN killcheck%=l 
NEXT J 
NEXT I 
END IF 

REM kill if no vacant or CO-filled A sites left 
IF Nrxn! (KK)=0 AND killcheck%=0 THEN GOTO labelSTART 
REM ****** SURFACE STILL ALIVE IF GET TO HERE **************** 

NEXT KK 

REM If allowed to complete FOR-NEXT loop, KK becomes M+l, so reset it to 
REM equal M. 

KK=KK-1 

REM average rates calculated after startup period 

Nrxn!(M+l)=0 'place holder for the ave of (KK+l-Avestart ) Nrxn! 

FOR KS= (startup%-rl) TO KK 

Nrxn! (M+l)=Nrxn! (M+1)+Nrxn! (KS ) / (KK-startup% ) 

NEXT KS s 
LOCATE 22 
PRINT fileout? 

PRINT "Ave . Rate =", Nrxn ! (K+l ) 

REM ***** check to see if this is the highest rate so far ******** 
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IF Nrxn! <M+1) >maxrace! THEN 
max rate ! =Nrxn ! (M+l) 

FOR 1=1 TO Ispan% 

FOR J=1 TO Jspan% 

maxsite% (I,J)=site% (I, J) 

NEXT J 
NEXT I 
BREAK OFF 
GOSUB suboutput 
END IF 

GOTO labelSTART 
labelBREAK: 

BREAK OFF 
GOSUB submaxgraph 
LOCATE 22 
labelASK : 

PRINT "Continue Program? <y or Y OR n or N): n 
INPUT MY $ 

IF MY$="y" OR KY$="Y" THEN GOTO labelSTART 

IF MY$<>"n" AND KY$<>"N" THEN CLS : LOCATE 12: GOTO labelASK 

OPEN f ileout $ FOR APPEND AS #1 

PRINT #l,"Total runs started in this set = ",runnum! 

PRINT #1,"" 

PRINT #1 , "END C? DATA" 

CLOSE #1 

END 

subranO : 

REM This routine is similar "RANO" in the book "Numerical Recipes (FORTRAN) 
REM W. H. Press, et ai. " 

REM It "scrambles" the "AND" numbers to cut down on pair correlations. 

REM I have also adeem a random re-seeding of "RND" to increase its repetitio 
REM period 

REM from originally cn the order of 5EQ7 calls (11,000 MC steps for 1092 REM 
REM array) . 

IF IFF=0 THEN 
IFF=1 

REM re-seed "AND" randomly with random seed to to increase its 
REM repetition period 
ZR%=32767*?ND 
IF RND< . 5 THEN ZR%=-ZR% 

RANDOMIZE ZR% 

REM the next lines "set up" the scrambling array for RANO 
FOR JR=i TO 97 
dum! =?irD 
NEXT JR 
FOR JR=1 TO 97 
VR! (JA)=RND 
NEXT JR 
YR ! =RND 

END IF 

JR=1+INT (97 ! * YR ! ) 

YR!=VR! (JR) 

RANO ! =YR ! 

VR! (JR) =RND 
RETURN 
subrandom: 

REM initialize for empty random A-B site distribution 
labelSRl : 
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NA% = 0 

FBsite ! =1- (NAbase%/ (l0span%*J0span%) ) 

FOR 1=1 TO I C span! 

FOR J=1 TO J0span% 

G0SU3 sub ran 0 

IF RANG ! >FBsite ! THEN 

siteO% (I, J) =0 :NA%=NA%+1 'site is empty A site 
GOTO labelSRO 
END IF 

sited (I, J) =10 T site is empty B site 

labelSRO: 

NEXT J 
NEXT I 

IF NA%ONAbasel THEN GOTO labelSRl * fixes no. A sites in base array 
RETURN 
suboldf ile : 

CLS 

labelASKl : 

filename$=FILES$ (1 , "TEXT" ) 

IF filename$= ,,n GOTO labelASKl 
PRINT "Input File is filenames 
OPEN filename$ FOR INPUT AS #1 
PRINT Mn 
MX=2 5 

FOR 1=1 TO MX 

LINE INPUT # 1 , dummy $ 

NEXT I 

LINE INPUT #1, na>:rate$ 

PRINT maxrate? 

MX=5 

FOR 1=1 TO MX 

LINE INPUT # 1 , dummy $ 

PRINT dummy $ 

NEXT I 

PRINT "hit any key to continue" 

INPUT MZZ 

FOR 1 = 1 TO I Ospar.% 

FOR J=1 TO JO span! 

INPUT #l,site0% (I, J) 

IF sitevd C, J) =1 THEN site0% (I, J)=0 
IF siteO* ( I , J) =2 THEN siteO % ( I , J) =0 
IF siteC% C, J) =12 THEN siteO % ( I , J) =1 0 
PRINT "slued (I, J) =", sited (I, J) 

NEXT J 
NEXT I 
CLOSE #1 

PRINT "enter any key to continue" 

INPUT MX 
CLS 

FOR IM=0 TO ( expander* -1 ) 

FOR JM=0 TO (expander* -1) 

FOR 1=1 TO I0span% 

FOR J=1 TO JOspanl 

REM set coordinates for "ovals" (circles here) in PAINTOVAL calls 
rectangle* (0) = { (JM* JOspanl) +J-1) *12 
rectangle* (1) = ( <IM*l0span% ) +1-1) *12 
rect angle % (2) - ( {JM* J0span% ) + J) *12 
rectangle% (3) = ( (IM*l0span%) +1) *12 
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dot% (0)=( ( ( JM* I Ospan! ) + J-l ) *12) +5 
dot% (1) = { { (IM*I0span%) +1-1) *12) +5 
dot % (2) = ( { (JM*J0span%) +J) *12) -5 
dot % (3) = ( { (IM*lOspan%) +1) *12) -5 
REM show current maximum rate site configuration 
IF siteO% (I, J) =0 THEN ! site is empty A site 
CALL PEN? AT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL {VARPTR ( rectangle% (0) ) ) 

GOTO labelOGRl 
END IF 

IF siteC% (I, J) =1 THEN ’site is CO ads on A site 
CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR {rectangle! (0) ) ) 

GOTO labelOGRl 
END IF 

IF siteO! (I, J) =2 THEN ’site is O ads on A site 
CALL PENP AT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (rectangle! (0) ) ) 

GOTO labelOGRl 
END IF 

IF siteO! (I, J) =10 THEN ’site is empty B site 
CALL SRASEOVAL (VARPTR (rectangle! (0) ) ) 

CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (dot % (0) ) ) 

GOTO labelOGRl 
END IF 

IF siteO! (I, J) =12 THEN ’site is O ads on B site 
CALL ERASEOVAL (VARPTR (rectangle! (0) ) ) 

CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (dot! (0) ) ) 

GOTO labelOGRl 
END IF 

PRINT "error: in subgraphics with unknown site configuration" 

REM shouldn't reach this point 
labelOGRl : 

NEXT J 
NEXT I 
NEXT JM 
NEXT IM 
LOCATE 18 
PRINT maxrate$ 

PRINT "" 

PRINT "hit any key zo continue: 

INPUT MX 
RETURN 
submanual : 

labelMAN : 

CLS 

PRINT "This is for base array =”, I Ospan! , JOspan! 

PRINT "With # A sites =",NAbase! 

PRINT "Input rows down first column, then second, then third, " 

PRINT "A's are 0, B’s are 10:" 

basecheck!=0 

FOR 1=1 TO I Ospan! 

FOR J=1 TO JOspan! 

INPUT siteO! (I, J) 

IF sited (I, J) =0 THEN basecheck!=basecheck!+l 
NEXT J 
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NEXT I 
CLS 

FOR IM=0 TO ( expander % -1 ) 

FOR JM=0 TO ( expander %-l ) 

FOR 1=1 TO I Ospan% 

FOR J=1 TO JOspan! 

REM set coordinates for "ovals" (circles here) in PAINTOVAL calls 
rectangle! (0) = ( ( JM* JO span!) +J-1) *12 
rectangle! (1) = ( (IM* IOspan!) +1-1) *12 
rectangle! (2 ) = ( (JM* JO span! ) + J) *12 
rectangle! (3) = ( (IM* iOspan! ) +1 ) *12 
dot! ( 0 ) = ( ( (JM*lOspan!)+J-l) *12) +5 
dot! (1) = ( ( (IM* IOspan!) +1-1) *12) +5 
dot! (2) = ( ( (JM* JO span!) +J) *12) -5 
dot! ( 3) = ( ( (IM* 10 span!) +1) *12) -5 
REM show current maximum rate site configuration 
IF siteO! (I, J) =0 THEN ‘site is empty A site 
CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR ( rectangle! (0) ) ) 

END IF 

IF siteO! (I, J) =1 0 THEN 'site is empty B site 
CALL ERASEOVAL (VARPTR (rectangle! (0) ) ) 

CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (dot! (0) ) ) 

END IF 
NEXT J 
NEXT I 
NEXT JM 
NEXT IM 
LOCATE 18 

IF NAbase!<>basecheck! THEN PRINT "NOT THE CORRECT NUMBER OF A SITES' 
PRINT "Enter 2 to revise, any other to accept:" 

INPUT MX 

IF MX=2 GOTO iabelMAN 
RETURN 
subexpand: 

NAsite ! =expander! * expander !*NAbase% 

I span%=l0 span %* expander! 

Jspan %= JO span ! * expander % 

FOR IM=0 TO (expander!-!) 

FOR JM=0 TO (expander! -1 ) 

FOR 1=1 TO IOspan! 

FOR J=1 TO J0span% 

site! ( CM* IOspan!) +1, ( JM* JOspan! ) + J) =site0% (I, J) 

NEXT J 
NEXT I 
NEXT JM 
NEXT IM 
RETURN 
suboutput : 

OPEN f ileout$ FOR OUTPUT AS #1 
PRINT #1 , f ileout $ 

PRINT #1, "" 

PRINT' #1 , "Runs tried to this point =",runnum! 

PRINT #1,"" 

PRINT #1 , "Initial state (l=random) = ", IX 
PRINT #1, "" 

PRINT # 1 , "NN (0 means diagonal actions not allowed) = " , NN 
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PRINT #1,"" 

PRINT #1 , "B3 (0 for 02 ads on AB pairs, 1=AB & BB, 2=AB & AA, 3=any" 
PRINT #1, "pair, 4=BB pairs) = ”,BB 
PRINT #1,”" 

PRINT #1, "kcoadsorb = ",kcoadsorb! 

PRINT #1 , "kcodesorb = ",kcodesorb! 

PRINT #1, "koxadsorb = ",koxadsorb! 

PRINT #1 , "kreact ! = ",kreact! 

PRINT #1, "yco!=",yco! 

PRINT #1,"” 

PRINT #1, "The total number of sites =",NStotal! 

PRINT #l,"The number of A sites =",NAsite! 

PRINT #1, "FAactual ! = ",FAactual! 

PRINT #1, "FBactual ! (compare to FBsite! for IX = 1) = ",FBactual! 
PRINT #1, "krxni =", krxn ! 

PRINT #1,"*’ 

PRINT #l,"Maxrate! is no. reactions per CO collision with an A site" 
PRINT #1, "averaged over last 24 of 28 X 100 MC steps" 

PRINT #1,"" 

PRINT #1, "Maxrare! ",maxrate! 

PRINT #1, "" 

PRINT #1, "The following is the site configuration, site0% (I, J) , " 
PRINT #1," listed as all JOspanl J's for 1=1, then all J0span% J's" 
PRINT #1," for 1=2, etc., up to I=I0span%" 

PRINT #1, "IOspar.% = ", l0span% 

PRINT #1 , " JOspar.l =",J0span% 

FOR 1=1 TO iOspanl 

FOR J=1 TO JOspanl 

PRINT el, siteO% (I, J) 

NEXT J 
NEXT I 
PRINT #1,"" 

CLOSE #1 
RETURN 
subpairpick : 

REM this picks a random pair of sites 
REM pick randcm site, site%(x,y) 

GOSUB subranO 

x=INT ( (Ispan%*RANO ! ) +1 ! ) 

GOSUB subranO 

y=INT ( (Jspar.%*RANO ! ) +1 ! ) 

REM pick neighboring site 

IF NN=0 THEN GOTO labelDD 

REM the following 7 lines are for: diagonal neighbor actions are OK 

labelNN: 

GOSUB subranO 
xn%=x+INT (3* RANO ! ) -1 
GOSUB subranO 
yn!=y+INT (3*RAN0 ! ) -1 

IF xn%=x AND yn%=y THEN GOTO labelNN 'neighbor can't be itself 
GOTO labe!599 

labelDD : 

REM the following 18 lines are for: diagonal actions not allowed 

GOSUB subranO 

IF RANO ! < . 5 THEN GOTO label99 
xn%=x 

GOSUB subranO 
IF RANO! <.5 THEN 
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yn%=ytl 

ELSE 

yn%=y-l 
END IF 

GOTO labe!999 

labe!99 : 

yn%=y 

GOSUB subranO 
IF RANO ! < . 5 THEN 
xnl=x-ri 

ELSE 

xn%=x-l 
END IF 
label999 : 

REM Set indices correctly for neighbors over array boundary. 

REM Corresponds to the surface as a torus or as a semi-infinite flat 
REM surface formed of a periodic array of the sub-array computed here. 
IF xn%=0 THEN xn%=Ispan% 

IF xn%>Ispar% THEN xn%=l 
IF yn%=0 THEN yn%= Jspan% 

IF yn%>Jspanl THEN yn%=l 

RETURN 
subcoadsorb : 

GOSUB subranO 

IF RANO ! <kcoadsorb! THEN site%(x,y)=l 
RETURN 
subcodesorb: 

GOSUB subranO 

IF RANO ! <kcodesort ! THEN site%(x,y)=0 
RETURN 
suboxadsorb : 

GOSUB subranO 

IF RANO ! <koxadsorb ! THEN 

site! (x, y) =site* y) + 2 
site% (xn%, yn% ) =site% (xn%, ynl ) +2 
END IF 
RETURN 
subrxn : 

GOSUB subranO 

IF RANO !>=. 4052345? THEN GOTO labelRXN 
Nrxn ! (KK)=Nrxn! (XK)+1! 

IF site% (x, y) =1 THEN sire%(x,y)=0 
IF site%(x,y}=2 THEN sire% (x , y ) =0 
IF site% (x, y) =12 THEN si te% (x, y ) =1 0 
IF site% (xn*,yr*)=l THEN site% <xn% , yn% ) =0 
IF site% (xnl,yn%)=2 THEN site! (xn%, yn% ) =0 
IF site% (xn%, yr%) =12 THEN site% (xn% , yn%) -1 0 
labelRXN: 

RETURN 
subgraphics : 

FOR 1=1 TO Ispanl 
FOR J-l TO Jspan% 

REM set coordinates for "ovals" (circles here) in PAINTOVAL calls 
rectangle* (0) = (J-l ) * 12 
rec tangled (1)=(I-1)*12 
rectangle* (2) =J*12 
rect angle % (3) =1*12 
dot% (0) = ( (J-l) *12)^5 
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dot% (1 ) = ( (1-1) *12 ) +5 
dot% (2) = ( J*12) -5 
dot % (3) - (1-12 ) -5 
REM update site graphics 

IF site%(I,J)=0 THEN 'site is empty A site 
CALL ERASEOVAL (VARPTR ( rectangle% (0) ) ) 

GOTO labelGRl 
END IF 

IF site%(I,J)=l THEN ‘site is CO ads on A site 
CALL PEN? AT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR ( rectangle% (0) ) ) 

GOTO labelGRl 
END IF 

IF site!(I,J)=2 THEN 'site is O ads on A site 
CALL PENPAT (VARPTR (shade! (0) ) ) 

CALL PAINTOVAL (VARPTR (rectangle! (0) ) ) 

GOTO labelGRl 
END IF 

IF site! ( I , J) =1 0 THEN 'site is empty B site 
CALL ERASEOVAL (VARPTR ( rectangle! (0) ) ) 

CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (dot! (0) ) ) 

GOTO labelGRl 
END IF 

IF site! (I , J) =12 THEN 'site is O ads on B site 
CALL PENPAT (VARPTR ( shade! (0) ) ) 

CALL PAINTOVAL (VARPTR (rectangle! (0) ) ) 

CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (dot! (0) ) ) 

GOTO labelGRl 
END IF 

PRINT "error: in subgraphics with unknown site configuration" 

REM shouldn't reach this point 
labelGRl : 

NEXT J 
NEXT I 
RETURN 
submaxgraph : 

CLS 

FOR 1=1 TO Ispar.% 

FOR J=1 TO Jspan! 

REM set coordinates for "ovals" (circles here) in PAINTOVAL calls 
rectangle! (0) = (J-l) *12 
rectangle! (1) = (1-1 ) *12 
rectangle! (2)=J*12 
rectangle! { 3 ) =1*12 
dot%(0)=( (J-l)*12)+5 
dot! (1) = { (I-l)*12)+5 
dot! (2) = { J*12) -5 
dot! (3) = (1*12) -5 

REM show current maximum rate site configuration 
IF maxsite! ( I , J) =0 THEN 'site is empty A site 
CALL PENPAT (VARPTR (black! (0) ) ) 

CALL PAINTOVAL (VARPTR (rectangle! (0) ) ) 

GOTO labelMGRl 
END IF 

IF maxsite! (I, J) =1 THEN 'site is CO ads on A site 
CALL PENPAT (VARPTR (black! (0) ) ) 
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CALL PAINTOVAL (VARPTR (rectangle% (0) ) ) 

GOTO labelMGRl 
END IF 

IF maxsite% (I, J) =2 THEN 'site is O ads on A site 
CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR ( rectangle% ( 0 ) ) ) 

GOTO labelMGRl 
END IF 

IF maxsite* ( I , J) =1 0 THEN 'site is empty B site 
CALL EPASEOVAL (VARPTR (rectangle% (0) ) ) 

CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR (dot% (0) ) ) 

GOTO labelMGRl 
END IF 

IF maxsite* (I, J) =12 THEN 'site is O ads on B site 
CALL EPASEOVAL (VARPTR (rectangle% (0) ) ) 

CALL PENPAT (VARPTR (black% (0) ) ) 

CALL PAINTOVAL (VARPTR (dot% (0) ) ) 

GOTO labelMGRl 
END IF 

PRINT "error: in subgraphics with unknown site configuration" 

REM shouldn't reach this point 
labelMGRl : 

NEXT J 
NEXT I 
LOCATE 18 

PRINT filecutS 

PRINT "Total number runs started = ", runnum! 

PRINT "Current cptimal surface (big black dots are A sites, small dots" 
PRINT "are B sices)" 

PRINT "Max. Race =",maxrate! 

RETURN 
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